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technology  can  be  considered  commercially 
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decisions  and  control  operations  that  need  to  be  made  only  once  for  each  data  packet.  Thus,  data 
throughput  through  a  single  router  can  be  increased  over  the  raw  electronic  processing  speed  by  a  factor 
equal  to  the  number  of  bits  in  the  packet,  i.e.,  several  orders  of  magnitude.  We  expect  our  concept  to 
provide  a  user  throughput  in  excess  of  30  Gb/s  and  an  overall  network  throughput  in  excess  of 
1012  bits/s. 

During  the  past  three  years,  SRI  International  investigated  many  issues  related  to  the  practical 
development  to  an  optical  packet  switch  using  rare-earth  doped  materials.  Specifically  SRI  focused  on 
three  major  research  areas.  (1)  Investigating  new  rare-earth  storage  materials  that  can  be  used  directly  at 
communication  wavelengths,  specifically  at  1.5  pm.  (2)  Demonstrating  optical  packet  buffering  at  this 
wavelength,  building  a  compact  demonstration  unit  using  optical  fiber  input  and  output.  (3)  Examining 
the  use  of  time-domain  holography  for  massively  parallel  optical  data  buffering/storage  and  determining 
the  data  input/output  transfer  rate  and  BER  of  the  system. 

Er3+:YAG  and  Er3+:Y2SiC>5  have  been  identified  as  suitable  candidates  for  data  buffering  at 
1.5  pm.  The  dephasing  time  has  been  measured  as  50  ns  for  Er3+:YAG  and  120  ps  for  ErJ+:Y2SiOs, 
with  the  latter  measurement  obtained  under  a  magnetic  field  of  6G.  The  material-limited  time-bandwidth 
products  for  both  materials  are  around  120,000,  which  suggests  a  potential  storage  density  of 
120  bit/focal  spot.  The  storage  time  for  three-pulse  echoes  has  also  been  measured  for  the  materials  and 
is  between  1  and  30  mn,  depending  upon  the  laser  power  used  in  recording  and  retrieval. 

In  a  separate  experiment  using  a  0.8-mm  spot  size  (FWHM=0.8  mm),  we  successfully  stored 
1020  bits  of  data  (306  ones  and  714  zeros)  in  a  40  MHz  channel  spectral  width.  The  total  length  of  the 
data  was  51  ps,  approximately  2/5  of  the  material's  phase  memory  time.  The  retrieved  data  intensity, 
therefore,  is  expected  to  decay  substantially  (approximately  one  order  of  magnitude)  over  the  length  of 
the  string.  To  prevent  this  fast  decay,  we  swept  the  intensity  of  the  input  data  exponentially  upwards  so 
that  the  entire  data  string  can  be  resolved  clearly  on  a  single  intensity  scale  with  a  good  signal-to-noise 
ratio.  Data  recall  can  take  place  at  any  time  after  storage.  This  attribute  made  this  technique  also  useful 
as  a  programmable  optical  delay  line. 

SRI  has  also  demonstrated  page-formatted  data  storage  using  the  time -domain  approach.  One 
hundred  pages  of  data,  each  of  which  contained  3,360  bits,  have  been  recording  in  Er3+:Y2Si05  at  a 
single  spatial  location  by  using  wavelength  multiplexing.  The  digital  data  were  recorded  and  accurately 
retrieved  at  a  peak  rate  of  300  Mbits  per  second  without  the  use  of  an  error  correcting  code.  The 
system's  bit  error  rate  is  about  10*7. 

The  experimental  results  from  this  project  clearly  demonstrate  the  practical  feasibility  of  using 
hole  burning  materials  for  optical  packet  switching  using  the  stimulated  echo  approach. 
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FIGURES 


1 .  Optical  random  access  memory  (ORAM) .  5 

2.  Optical  routing  switch,  (a)  General  schematic;  (b)  source,  read,  write,  and 

destination  packet  directions .  6 

3.  Energy  level  diagram  of  Er3+:YAF  for  the  4Ii5/2  and  4Ii3/2  state .  8 


4 .  (a)  Two-pulse  echo  intensity  in  Er3+: YAG  at  -7  K  under  an  external  magnetic 
field  of  100  G  as  a  function  of  change  of  excitation  pulse  separation.  Open 
circles  are  data  and  the  solid  line  is  a  simple  exponential  fit.  (b)  Stimulated 
photon  echo  intensity  in  Er3+:YAG  as  a  function  of  the  separation  between  the 
second  and  the  third  excitation  pulses.  The  results  are  obtained  at  ~1  K  under 

an  external  magnetic  field  of  3000  G .  10 

5 .  Storage  of  1020  bits  of  data.  The  string  is  the  ASCII  representation  of  "Time- 
Domain  Holography,  Dynamic  Optical  Memory,  Reprogrammable  Optical 
Delay  Lines  at  1.5  microns"  encoded  with  the  (1.7)  code.  An  exponential 
intensity  ramp  is  used  to  compensate  to  the  large  decrease  of  the  intensity  of 
the  retrieved  data  over  50  (is.  The  peak  power  is  measured  to  be  -0.25  pW. 

All  of  the  data  traces  shown  in  this  figure  represent  single  recording  and 
retrieval  events  and  are  detected  by  an  avalanche  photodiode  module  with  a 

gain  of  9.4  V/pW  at  1.55  pm .  13 

6 .  (a)  Input  data  image  after  transmitted  through  the  Eu3+:  Y2Si05  sample. 

(b)  Reconstructed  image  after  being  stored  in  Eu3+:Y2SiC>5  for  200  ps. 

(c)  Image  reconstructed  under  the  same  condition  as  that  of  (b)  but  with  weak 
overall  intensity,  (d)  Image  reconstructed  5  minutes  after  the  completion  of 
recording.  In  (a)-(c),  the  data  pulse  length  was  1 1.6  ps,  and  that  in  (d)  was  50 


ps.  The  former  corresponds  to  an  I/O  rate  of  300  megabits  per  second  while 

the  latter  corresponds  to  67  megabits  per  second .  16 

7.  (a)  Histogram  showing  the  intensity  distributions  of  0’s  and  l’s  that  were 
retrieved  from  the  time-domain  memory,  (b)  Normalized  histogram  obtained 

after  applying  the  local  threshold  technique .  17 

8 .  (a)  Digitized  photograph  displayed  using  its  original  binary  data,  (b)  Same 
photograph  obtained  after  the  binary  data  being  transferred  to  and  then  recalled 

back  from  the  time-domain  memory .  20 
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There  is  a  growing  demand  for  high  throughput  data  networks  in  both  military  and 
commercial  applications,  such  as  distribution  of  intelligence  information  (data,  images)  for 
the  military  and  communication  of  multimedia  services  (voice,  video,  image,  and  data)  for  the 
commercial  sector.  Low-cost  desktop  workstations  with  more  than  100  MIPS  of  processing 
power  are  currently  being  used,  and  workstations  with  -1000  MIPS  will  be  available  in  a  few 
years.  The  image-based  interaction  and  collaboration  applications  running  on  these  workstations 
will  impose  enormous  demands  on  network  capacity.  For  example,  transmission  of  a  single 
high  density  television  (HDTV)  channel  requires  a  data  rate  in  excess  of  40  Mbits/s,  even  with 
compression.  Networks  of  the  next  decade  will  need  to  support  hundreds  or  even  thousands  of 
such  users.  It  is  clear  that  the  gigabit  network  architectures  being  developed  today  will  not  be  able 
to  meet  the  capacity  requirements  in  the  next  decade. 

With  fiber  optics  as  the  physical-layer  technology,  the  limiting  factor  on  the  bandwidth 
of  a  network  is  not  the  transmission  medium  but  the  electronic  and  photonic  devices  at  the  end¬ 
points.  A  single-mode  optical  fiber  operating  at  a  wavelength  in  the  neighborhood  of  1550  nm 
has  a  capacity  of  approximately  25,000  GHz,  while  electronic  processing  speeds  are  limited  to 
a  few  Gbits/s.  This  speed  mismatch  can  be  expected  to  prevail  for  a  long  time.  This  limitation 
led  researchers  to  investigate  all-optical  networks,  where  the  path  between  two  communicating 
end-points  does  not  involve  conversion  to  electronic  signals  for  processing. 

If  optical  networks  are  to  be  used  widely,  packet  switching  will  need  to  be  supported,  as 
well  as  methods  for  interconnection  of  these  networks.  Packet  switching  allows  more  efficient 
utilization  of  network  capacity  than  circuit  switching  because  of  the  statistical  multiplexing  of 
bursty  users.  More  importantly,  packet  switching  allows  multiple  upper-level  protocols  and 
services  to  share  the  same  physical  network  paths  concurrently.  Such  support  for  multiple  con¬ 
current  traffic  streams  is  an  absolute  requirement  in  interactive  computing.  In  addition,  a  large 
number  of  local  area  and  wide  area  networks  in  the  military  and  commercial  sectors  will  be  packet 
switching  networks  to  which  future  optical  networks  must  be  connected.  Hence,  from 
the  Department  of  Defense's  communication  infrastructure  standpoint,  any  practical  long-term 
solution  to  the  utilization  of  optical  devices  for  networking  must  accommodate  packet  switching. 

To  illustrate  the  concept  of  optical  switching  using  the  stimulated  echo  concept,  we  use  a 
simple  diagram  to  describe  the  storage  and  retrieval  of  time  domain  data  and  then  describe  the 
switching  aspects  of  this  technique. 
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As  shown  in  Figure  la,  information  is  written  as  phase/amplitude  or  frequency  modu¬ 
lations  by  encoding  the  Data  and  Write  laser  beams,  which  cross  inside  the  sample,  at  times  tD 
and  tw,  respectively.  The  relative  directions  of  the  laser  pulses  are  shown  in  Figure  lb.  The 
overlap  of  these  two  beams  causes  the  dopant  materials,  such  as  rare-earth  ions  doped  in  oxide 
crystals,  to  Fourier  transform  the  time-domain  data  and  store  the  information  as  a  frequency 
grating  at  low  temperatures.  Information  is  retrieved  sometime  later  by  illuminating  the  crystal  by 
a  single  Read  pulse  (which  also  may  be  modulated)  at  time  tR.  The  crystal  then  emits  a  laser-like 
coherent  beam  (the  Echo  beam)  that  is  a  replica  of  the  original  Data  beam  at  time  te  =  tR  +  (tD-tw)- 
This  random  access  memory  can  be  dynamically  configured  with  word  sizes  ranging  from  one  bit 
up  to  106  bits.  Memory  time  of  up  to  a  day,  memory  density  of  1010  bits/cm3,  and  single 
pixel/channel  I/O  throughput  in  excess  of  100  Mbits/s  have  been  experimentally  demonstrated. 
Figures  lc  through  le  show  the  experimental  storage  and  retrieval  of  0.4  kbits  of  information  in  a 
single  spatial  spot  in  Eu^SiOs.  We  have  recently  demonstrated  storage  of  2.0  kbits  of  data  in  a 
single  spatial  spot  in  Eu:Y2Si05. 

ANGULAR  ROUTING,  SWITCHING,  AND  MULTICASTING 

The  storage  medium  stores  not  only  the  information  carried  by  the  Data  and  Write  beams 
w  ■■■■)  — )  ■ )  ■ 

but  also  the  relative  angle  of  the  beams.  It  can  be  shown  that  kecho  =  +  kR  -  kw,  where 

"kecho.  kw»  kR.  and  kD  are  the  propagation  wave  vectors  of  the  retrieved  data  (Echo),  Write, 

Read,  and  Data  pulse  train.  This  phase  matching  condition  is  exploited  to  produce  angular 

routing,  switching,  and  multicasting  operation.  Figure  2a  shows  a  schematic  outline  of  the  optical 

routing  switch.  Figure  2b  shows  how  buffered  packet  switching  using  the  phase  matching 

conditions  discussed  above  can  be  accomplished.  Assume  that  the  packet  exiting  an  input 

(source)  fiber  port  i  needs  to  be  space-switched  to  an  output  (destination)  port  j.  The  destination 

port  information  is  contained  in  the  header  information  of  the  packet.  This  header  information  is 

optically  decoded  as  discussed  below  and  the  output  port  determined.  The  Data  and  Write  pulses 

intersect  inside  the  sample  at  angle  Tt-0.  The  destination  port  j  is  located  at  an  angle  0  with 

respect  to  the  source.  The  input  packet  emerging  from  the  source  fiber  impinges  on  the  crystal 

after  passing  through  a  fiber  delay  line.  This  data  packet  is  first  stored  (buffered)  in  a  particular 

physical  location  in  the  crystal  by  addressing  the  location  first  by  a  Write  pulse  and  then  by  the 

data  packet.  This  procedure  stores  not  only  the  data  packet  information  but  also  the  destination 

angle  0.  The  stored  packet  may  be  routed  and  switched  out  of  the  memory  by  addressing  the 

memory  location  with  a  short  Read  pulse  that  travels  opposite  to  the  input  packet.  This  Read 

pulse  causes  the  memory  to  coherently  emit  (transmit)  a  packet  (echo)  that  mimics  the  input  data  at 

an  angle  0  with  respect  to  the  input  data  direction,  as  shown  Figure  2b.  The  multicasting 

operation  is  easily  accomplished  by  simultaneously  sending  two  or  more  Write  pulses  at  angles 

01,  02,  and  so  on,  during  the  recording.  The  packet  thus  stored  can  be  readily  switched  to  the 

predetermined  destination(s)  by  using  a  single  Read  pulse. 
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(a)  Temporal  schematic  of  first-in-first-out  (FIFO)  data  memory. 


(b )  One  of  the  many  possible  relative  directions  of  the  write,  data, 
read,  and  echo  pulses. 


(c)  Input  pulse  sequence  and  echo  signal.  Full  horizontal  scale  - 
102.4  jis.  The  duration  of  the  data  pulse  is  10  jis  (400  bits). 
The  input  pulses  are  attenuated  by  -KT4. 


(d)  Input  data  pulse  (P2)  on  an  expanded  time  scale.  Full  horizontal 
scale  -  2.4  |is. 


(e )  Retrieved  data  pulse  (echo)  on  an  expanded  time  scale. 


CAM-330581-167 


Figure  1 .  Optical  random  access  memory  (ORAM). 
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Figure  2.  Optical  routing  switch,  (a)  General  schematic;  (b)  source,  read, 
write,  and  destination  packet  directions. 
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RESULTS 


TIME  DOMAIN  MATERIALS  AT  1.5  |im 

Unlike  other  optical  memories,  TDOM  is  a  high-speed  signal  processor  by  nature. 
Temporal  data  strings  stored  at  different  times  can  be  correlated  or  convoluted  simultaneously 
with  an  incoming  data  string  in  the  memory  at  bandwidths  up  to  several  GHz  with  no  additional 
processing  time  required  [1-8].  This  rapid  in-memory  processing  feature  combined  with  the 
demonstrated  high  storage  capacity  makes  TDOM  also  a  superior  candidate  for  applications 
ranging  from  high-speed  optical  packet  switching,  context-based  retrieval,  matched  filtering,  data 
encryption,  to  pattern  recognition. 

In  spite  of  the  rapid  progress  in  the  development  of  optical  memories  and  all-optical 
processors,  to  our  knowledge  there  exists  no  model  directly  operating  at  the  optical 
communication  wavelengths.  The  current  optical  devices  require  electrical  to  optical  signal 
conversions  which  certainly  reduce  the  data  communication  speed.  To  meet  the  rapidly  increasing 
demand  for  high  speed  and  high  capacity  optical  communications,  it  becomes  necessary  to 
eliminate  such  electrical-optical  conversions.  However,  it  remains  a  challenge  to  find  a  suitable 
medium  to  store  and  process  optical  signals  at  the  infrared  wavelengths  of  fiber  communication. 

SRI  has  explored  the  feasibility  of  developing  a  TDOM  that  can  be  operated  directly  at  the 
communication  wavelengths.  We  choose  the  4Ii5/2  -  4Ii3/2  transition  in  Er3+:YAG  as  a  candidate 
for  the  present  study  and  examine  some  of  its  spectroscopic  properties  relevant  to  time-domain 
applications. 

The  erbium  ion  has  an  odd  number  of  f-electrons,  giving  rise  to  Kramers  degeneracy, 
which  can  only  be  lifted  by  external  magnetic  field.  We  resolved  six  lines  associated  with 
transitions  from  the  lowest  Stark  level  in  the  41 15/2  state  to  six  Stark  levels  of  the  4Ii3/2  state  (see 
Figure  3).  Their  wave  numbers  in  order  are  6549,  6599,  6606,  6786,  6806,  and  6886  cm-1.  The 
absorption  line  at  6527  cm4  at  ~7  K  is  a  hotband  from  the  second  lowest  Stark  level  (~22  cm*1)  in 
the  4Ii5/2  state  to  the  lowest  level  of  the  4Ij3/2  state.  At  higher  temperature,  we  also  observed 
hotband  from  higher  Stark  levels  of  the  4Ii3/2  state.  The  energy  level  diagram  of  4Ii5/2  and  4Ii3/2 
corresponding  to  our  data  is  shown  in  Figure  3.  Our  results  are  in  general  consistent  with  earlier 
work  on  Er3+:YAG  by  Ashurov  and  coworkers  [9],  who  also  observed  only  six  lines  for  the 
4Il5/2 ' 4Ii3/2  transition  at  4  K.  The  data  reported  by  Koningstein  and  Geusic  [10]  indicates  an 
additional  absorption  line  ~91  cm*1  above  the  highest  line  in  our  data,  which  may  be  too  weak  for 
us  to  detect. 
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Figure  3.  Energy  level  diagram  of  Er^iYAF  for  the  4Ii 5/2  and  4Ii3/2  state. 
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The  inhomogeneous  broadening  of  the  transition  was  estimated  to  be  -20  GHz  from  the 
measured  transmission  spectrum.  The  homogeneous  line  width  of  this  transition  was  measured 
by  recording  two-pulse  echo  intensity  as  a  function  of  excitation  pulse  separation.  A  simple 
exponential  fit  to  the  data  gives  a  homogeneous  dephasing  time  T2  of  -50  ns.  This  dephasing 
time  was  found  to  vary  little  between  1  and  7  K  under  a  magnetic  field  below  100  G. 

In  our  two-pulse  echo  experiment,  we  used  two  dye  lasers  with  10  GHz  spectral  width. 
The  peak  intensity  of  the  laser  inside  the  crystal  is  about  1  MW/cm2.  We  use  the  fluorescence 
lifetime  of  the  4Ii3/2  excited  state  to  calculate  the  laser  pulse  area  [11]  and  take  into  account  that 
the  laser  pulse  is  non-Fourier-transform-limited.  The  estimated  pulse  area  is  between  71/8  and 
7t/4.  This  area  of  <  1  was  also  supported  by  the  experimental  observation  that  the  echo  intensity 
was  linear  to  the  first  laser  pulse’s  intensity  and  quadratic  to  the  second.  We  also  found  that 
reducing  the  laser  power  by  a  factor  of  2  had  no  effect  on  the  decay  rate  of  the  echo. 

The  data  in  Figure  4a  were  taken  with  an  external  magnetic  field  of  -100  G  applied  along 
the  c-axis  of  the  crystal.  When  this  external  field  was  raised  from  zero  to  10  G,  the  echo  intensity 
quickly  increased  by  a  factor  of  three.  But  this  increase  became  saturated  when  the  field  was 
above  30  G.  On  the  other  hand  the  dephasing  time  did  not  change  upon  increasing  the  magnetic 
field  to  100  G.  When  a  magnetic  field  of  -3000  G  was  applied  to  the  crystal  (perpendicular  to  its 
c-axis),  we  observed  apparently  non-exponential  decay  of  echo,  but  the  dephasing  time  was  not 
lengthened  either. 

One  of  the  important  parameters  governing  the  performance  of  the  time-domain  memory  is 
the  highest  temperature  at  which  the  storage  sample  can  be  operated.  To  determine  this 
temperature,  we  measured  two-pulse  echo  intensity  at  different  sample  temperatures.  At  a  fixed 
excitation  pulse  separation  about  35  ns,  the  echo  intensity  gradually  decreased  by  only  a  factor  of 
two  from  7  to  1 1  K.  In  comparison  to  many  other  rare-earth  doped  crystals  except  Eu3+:Y2Si05 
at  its  7Fq-5D0  transition  [12],  this  temperature  dependence  of  Er3+:YAG  is  considerably  weaker. 

In  addition  to  the  two-pulse  echo  experiments,  we  also  examined  three-pulse  (or 
stimulated)  photon  echo  in  Er3+:YAG  since  it  is  closely  related  to  the  physical  processes  involved 
in  TDOM  and  related  optical  processors.  By  doing  so,  we  would  be  able  to  estimate  the  potential 
performance  of  Er3+:YAG  as  TDOM  at  1.527  pm,  such  as  the  maximum  storage  time,  which  is 
often  limited  by  spectral  diffusion  resulting  from  ion-host  and  ion-ion  interactions.  We  carried 
out  the  three-pulse  stimulate  echo  experiment  under  a  condition  very  similar  to  the  two-pulse  echo 
experiment.  We  fixed  the  separation  between  the  first  two  excitation  pulses  at  20  ns  and  varied 
that  between  the  second  and  the  third  pulses.  The  results  (shown  in  Figure  4b)  show  a 
complicated  decay  feature,  with  an  initial  fast  decay,  followed  by  a  slower  decay  and  then  another 
fast  decay.  At  a  delay  of  50  ps,  the  echo  intensity  decreased  by  a  factor  of  10  from  that  at  1  ps. 
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Figure  4.  (a)  Two-pulse  echo  intensity  in  Er^YAG  at  -7  K  under  an  external  magnetic  field  of  100  G  as 
a  function  of  change  of  excitation  pulse  separation.  Open  circles  are  data  and  the  solid  line  is  a  simple 
exponential  fit.  (b)  Stimulated  photon  echo  intensity  in  Er^YAG  as  a  function  of  the  separation  between 
the  second  and  the  third  excitation  pulses.  The  results  are  obtained  at  -1  K  under  an  external  magnetic 
field  of  3000  G. 
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This  decay  is  much  faster  than  the  fluorescence  lifetime  measured  earlier,  suggesting  the 
presence  of  rapid  spectral  diffusion  on  a  microsecond  time  scale.  One  possible  cause  of  this  fast 
diffusion  is  the  strong  Er-Er  dipolar  interactions  as  suggested  by  Macfarlane  et  al.  [13]. 

In  conclusion,  we  observed  for  the  first  time  two-pulse  and  three -pulse  echoes  at  1.527 
Jim  in  Er3+:YAG.  The  fast  dephasing  together  with  its  insensitivity  to  external  magnetic  field 
observed  in  our  experiments,  although  less  desirable,  pose  no  significant  challenge  to  the 
implementation  of  Er3+:YAG  as  a  time-domain  storage  medium  in  most  of  applications.  For 
example,  recently  developed  swept-carrier  technique  allows  the  storage  of  data  strings  that  are 
much  longer  than  the  dephasing  time  of  a  storage  material.  The  lifetime  of  the  spectral  grating 
measured  by  three-pulse  echoes,  although  short,  is  still  adequate  for  applications  such  as  dynamic 
optical  memory,  optical  packet  buffering,  and  optical  signal  processing.  The  present  work  opens 
up  the  possibility  of  developing  TDOM  and  TDOM-based  optical  processors  operated  directly  at 
the  communication  wavelengths.  Further  improvement  of  material  properties,  such  as  longer 
dephasing  time  and  less  spectral  diffusion,  may  be  accomplished  by  choosing,  e.g.,  different  host 
materials,  reducing  Er  concentrations,  and/or  applying  larger  external  magnetic  fields.  These 
efforts  have  shown  success  in  other  rare-earth  ions  such  as  Eu  and  Pr  ions.  However,  they  may 
not  be  applicable  to  Er  ions. 

ALL-OPTICAL  DYNAMIC  MEMORY  FOR  FIBER  COMMUNICATIONS 

In  time-domain  optical  memory,  temporal  information  is  stored  as  spectral  interference 
patterns  in  an  inhomogeneously  broadened  electronic  transition  of  rare-earth  ions  by  resonant 
excitations.  The  interference  patterns  (referred  to  as  spectral  gratings)  are  generated  by  using  a 
reference  pulse  known  as  the  write  pulse,  which  interferes  in  time -domain  or  equivalently  in 
spectral  domain,  with  the  data  sequence  to  be  stored.  Each  bit  of  information  is  recorded  as  an 
integral  part  of  the  spectral  gratings  which  extends  over  a  spectral  region  determined  by  the  data 
modulation  bandwidth.  These  gratings  in  a  carefully  chosen  material  can  persist  for  several  hours  to 
perhaps  several  days  at  cryogenic  temperatures  because  the  recorded  information  is  actually  stored 
as  population  modulation  of  hyperfme  levels  in  the  ground  state.  The  maximum  bandwidth  of  a 
data  string  is  limited  by  the  inhomogeneous  linewidth  of  an  electronic  transition  to  be  used,  which 
for  rare-earth  doped  solids,  is  typically  between  1.0  and  10  GHz. 

Temporal  information  recorded  in  the  memory  is  retrieved  by  illuminating  the  spectral 
gratings  with  a  read  pulse,  which  induces  coherent  radiation  known  as  photon  echoes.  In  order  for 
the  echo  to  mimic  the  temporal  data  string  originally  stored,  the  product  of  the  Fourier  transforms  of 
the  write  and  read  pulses  must  be  spectrally  flat  over  the  bandwidth  of  the  data  string.  This  feature 
of  rapidly  recording  optical  information  in  a  bit-serial  format  distinguishes  the  time-domain  memory 
from  other  optical  storage  techniques,  such  as  photorefractive  and  persistent  spectral  holebuming 
memories,  where  serial  data  must  be  first  converted  to  parallel  ones  before  recording  can  take  place. 
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Storage  of  1020  bits  of  data  is  shown  in  Figure  5.  The  storage  material  is  Er3+:YiSi05 
(0.005  at.  %  and  7  mm  thick).  Optical  information  is  stored  in  the  transition  of  the 

crystal,  which  has  two  optically  equivalent  sites,  one  at  1.5364  p.m  (site  1)  and  the  other  at  1.5388 
pm  (site  2).  Each  optical  site  has  an  inhomogeneous  broadening  of  -1  gigahertz.  To  increase  the 
length  of  a  data  string  that  can  be  stored  each  time  (which  is  determined  by  the  phase  memory  time 
of  the  transition),  we  place  the  sample  between  two  permanent  magnets  with  a  magnetic  field  (~  3 
kilogauss)  perpendicular  to  the  optical  beam  axis. 

A  number  of  issues  still  need  to  be  addressed  before  such  a  dynamic  memory  can  be  built. 
One  of  them  is  how  to  tailor  the  gain  of  the  optical  amplifier  so  that  the  retrieved  optical  waveforms 
are  identical  to  the  input  to  maintain  the  signal  fidelity  over  a  large  of  refreshing  cycles.  This  may  be 
accomplished  by  saturating  the  gain  medium  of  the  optical  amplifier,  for  example. 

MASSIVELY  PARALLEL  TIME-DOMAIN  OPTICAL  MEMORY 

Time -domain  holographic  memory  has  emerged  as  an  excellent  candidate  for  such  I/O 
intensive  applications.  This  memory  possesses  not  only  features  that  are  common  to  conventional 
spatial  holographic  recording,  such  as  high  capacity  and  high  degree  of  parallelism,  but  also  those 
that  are  unique  and  complementary  to  the  conventional  holography.  The  most  appealing  one  is  its 
fast  page  recording  speed  under  a  modest  laser  power.  Recent  experiments  [7]  show  that  an 
exposure  time  of  several  microseconds  and  a  peak  laser  power  of  less  than  200  mW  are  sufficient 
to  record  a  high-resolution  binary  image  with  the  time-domain  technique. 

The  spectral  interference  of  the  two  pulses  is  recorded  in  a  storage  medium  (which  are  often 
rare-earth  doped  solids  such  as  Eu3+:Y2Si05  or  Pr3+:  YAG)  by  laser  excitation.  Ions  such  as  Eu3+ 
and  Pr5+  that  are  responsible  for  data  recording  are  resonantly  excited  by  the  laser  pulses  to  form 
population  gratings  resembling  the  Fourier  spectra  of  the  pulse  pair.  These  gratings  in  a  carefully 
chosen  material  can  persist  for  several  hours  to  perhaps  several  days  at  cryogenic  temperatures 
because  the  recorded  information  is  actually  stored  as  population  modulation  of  hyperfine  levels  in 
the  ground  state  [7].  Since  the  grating  is  a  direct  result  of  resonant  excitation  and  requires  no 
secondary  process,  its  formation  is  instantaneous  and  requires  modest  laser  power,  the  data 
recording  is  therefore  inherently  fast  and  efficient,  giving  rise  to  high  date  transfer  rates,  fast  access 
times,  and  short  latency.  These  features  are  ideal  for  applications  as  a  dynamic  memory. 

The  storage  system.  The  material  chosen  to  demonstrate  this  novel  holographic 
recording  scheme  is  crystalline  Eu3+:Y2SiO5(0.1  at.  %  and  7  mm  thick).  We  use  its  Tq  -  5D0 
transition  at  Site  1  (-579.88  nm)  which  has  an  inhomogeneously  broadened  linewidth  of  -4  GHz 
and  a  dephasing  time  of  -2  ms.  The  product  of  the  two  parameters  gives  the  theoretical  upper 
bound  for  the  storage  capacity  per  spatial  location  as  a  result  of  the  presence  of  the  wavelength 
dimension.  This  material  has  two  optically  equivalent  sites  for  the  transition;  the  other  site  (Site  2) 
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Figure  5.  Storage  of  1020  bits  of  data.  The  string  is  the  ASCII  representation  of  “Time-Domain  Holography,  Dynamic 
Optical  Memory,  Reprogrammable  Optical  Delay  Lines  at  1.5  microns"  encoded  with  the  (1.7)  code.  An  exponential 
intensity  ramp  is  used  to  compensate  to  the  large  decrease  of  the  intensity  of  the  retrieved  data  over  50  ps.  The  peak 
power  is  measured  to  be  -0.25  pW.  All  of  the  data  traces  shown  in  this  figure  represent  single  recording  and  retrieval 
events  and  are  detected  by  an  avalanche  photodiode  module  with  a  gain  of  9.4  V/pW  at  1 .55  pm. 


has  a  comparable  linewidth  and  dephasing  time,  making  the  material  a  good  candidate  for  high 
capacity  applications.  During  the  experiments,  the  crystal  was  kept  at  a  temperature  of  about  4  K 
in  a  flowing  He-vapor  cryostat. 

Each  data  page  is  stored  in  a  distinct,  narrow  (-500  kHz)  frequency  channel  within  the 
inhomogeneously  broadened  absorption  line  by  wavelength  division  multiplexing  (WDM).  The 
separation  between  two  adjacent  channels  is  chosen  to  be  approximately  800  kHz.  WDM  has 
several  distinct  advantages.  Recording  and  retrieval  can  be  accomplished  without  use  of  any 
moving  parts  for  laser  beam  steering.  Therefore,  it  is  inherently  reliable  and  often  provides  short 
access  times.  In  addition,  WDM  is  expected  to  have  no  significant  effect  on  diffraction  efficiency. 
Each  data  page  can  be  recorded  in  an  identical  manner,  no  scheduled  recording  or  incremental 
recording  as  used  in  angular  division  multiplexing  [14]  are  needed,  increasing  considerably  the 
recording  speed. 

The  spatial  light  modulator  used  in  the  present  experiment  to  generate  two  dimensional 
digital  data  is  a  liquid  crystal  array  taken  from  a  commercial  projection  television.  The  array 
consists  of  480  x  440  pixels  and  has  a  low  fill  factor  of  -24%.  This  low  fill  factor  introduces  a 
large  insertion  loss  in  the  transmitted  beam.  We  measured  the  loss  for  the  zeroth-order 
transmission  to  be  -97%  at  580  nm. 

The  retrieved  data  was  captured  by  an  intensified  CCD  camera  which  has  an  array  of 
768  x  484  pixels.  The  captured  images  were  further  digitized  at  reduced  horizontal  resolution  of 
640  pixels  per  line  by  a  frame  grabber.  The  mismatch  between  the  CCD  and  SLM,  and  that 
between  the  frame  grabber  and  CCD,  as  well  as  the  low  spatial  resolution  of  the  intensifier  (-64 
lp/mm),  prevented  the  use  of  small  bit  dimensions  in  our  experiment.  We  chose  5x5  pixels  on 
the  SLM  to  represent  a  data  bit  and  used  a  total  of  3360  bits  (60  x  56  bits)  in  each  data  frame  to 
cover  a  center  area  equal  to  approximately  40%  (— 1.6  x  1.3  cm')  of  the  SLM.  The  data  frame  was 
imaged  (1:1)  on  to  the  CCD  array,  resulting  in  an  average  bit  size  of  10.7  x  8.4  pixels  after 
digitization.  We  used  these  mismatched  devices  because  they  were  readily  available  at  the  time  of 
the  experiments,  and  no  attempt  was  made  to  match  the  dimensions  of  the  devices.  Obviously, 
choosing  matched  SLM,  CCD  and  frame  grabber  would  allow  the  use  more  data  bits  per  frame, 
which  would  increase  the  storage  capacity  as  well  as  the  data  transfer  rate. 

The  write/read  pulses  used  in  the  present  experiment  were  10  pis  long  square  pulses  with  a 
peak  power  of  -200  mW  only.  These  pulses  were  biphase-encoded  with  the  5-bit  Barker  code 
[15]  to  increase  data  channel  width  to  -500  kHz.  The  reason  for  introducing  this  modulation  is  to 
make  our  system  more  tolerant  to  laser  wavelength  instability  and  hence  reduce  the  fluctuation  of 
echo  intensity.  The  length  of  the  data  pulse  used  in  our  experiment  varied  from  1 1.2  to  50  ps. 
This  data  pulse  was  also  biphase  modulated  by  using  pseudorandom  code  [16]  to  further  reduce 
the  echo  fluctuation  caused  by  laser  wavelength  instability. 

Both  the  SLM  and  CCD  operated  at  the  maximum  video  rate  (i.e.,  30  frames  per  second) 
and  the  resulting  digitized  data  was  stored  in  the  dynamic  memory  of  a  computer  at  the  same 
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speed.  Compared  to  1 1.6  (is  per  page  exposure  time  used  in  our  experiment,  this  video  rate 
obviously  limits  the  sustainable  I/O  speed  of  our  system  substantially.  This  limit  however  is  not 
fundamental  and  can  be  increased  to  match  the  speed  set  by  the  exposure  time. 

Analysis  of  page-formatted  binary  data.  In  digital  memories,  one  important 
parameter  determining  the  memory’s  performance  is  the  bit  error  rate  (BER),  which  measures 
how  often  an  error  bit  would  occur  on  average.  This  parameter  is  closely  related  to  other 
performance  parameters  of  the  system,  such  as  storage  capacity,  data  transfer  rates  and  data 
access  times.  For  a  given  system  in  general,  attempting  to  reduce  the  BER  often  leads  to  a 
decrease  of  storage  capacity  and/or  data  transfer  rate,  because  a  low  BER  requires  high  signal-to- 
noise  ratios.  For  parallel  recording,  it  means  the  need  of  a  large  modulation  transfer  function 
which  implies  the  use  of  large  bit  dimensions.  For  a  fixed  frame  size,  for  example,  it  translates 
into  a  decrease  of  the  data  transfer  rate. 

For  an  optical  storage  system  operating  at  high  I/O  speeds  with  low  laser  power 
illumination,  the  retrieved  data  are  often  at  a  signal  level  close  to  shot  noise  limited.  For  page- 
formatted  recording,  it  is  further  complicated  by  the  nonuniform  illumination  of  a  SLM.  Because 
of  these  undesired  conditions  associated  with  low  photon  budget,  such  a  system  often  experiences 
large  page-to-page  fluctuations  and  intra-page  variations. 

Figure  6a  shows  a  sample  input  page  used  in  data  recording.  Here,  spatial  variation  of 
intensity  is  apparent.  This  variation  reflects  the  spatial  beam  mode  of  the  addressing  laser, 
contributing  partially  to  the  spatial  signal  variations  within  a  page  discussed  earlier.  The  other 
contributions  include  shot  noise  resulting  from  the  quantum  nature  of  photons. 

The  page-to-page  fluctuation  is  often  more  noticeable  in  the  retrieved  data  because  it  is 
magnified  by  less  than  ideal  recording  and  readout  conditions.  Figures  6b  and  6c  show  two 
retrieved  pages  that  represent  the  brightest  and  dimmest  among  a  total  of  100  pages  recorded  in  an 
experiment.  They  are  single  readout  events  with  no  signal  averaging.  The  intensity  difference 
between  the  two  pages  is  >300%. 

Further  examining  the  histogram  of  all  100  pages  (totaling  336  kbits  of  data),  we  found 
that  the  intensity  distributions  of  the  l’s  (or  ON’s)  and  0’s  (or  OFF’s)  partially  overlap  as  shown 
in  Figure  7.  If  a  single  threshold  were  used  to  determine  all  the  data  bits,  it  would  lead  to  a  BER 
as  large  as  10'3,  not  suitable  for  any  storage  application. 

We  use  recently  determined  data  bits  to  obtain  a  local  threshold  for  the  next  bit  under 
consideration.  The  local  region  is  chosen  in  such  a  way  it  occupies  an  area  that  is  small  enough 
so  that  the  illumination  is  approximately  uniform,  yet  large  enough  to  have  a  sufficient  number  of 
bits  for  an  adequate  average.  Under  such  a  condition,  a  threshold  can  be  derived  from  the 


13 


(d) 

Figure  6.  (a)  Input  data  image  after  transmitted  through  the  Eu3+:Y2SiOs  sample,  (b)  Reconstructed 
image  after  being  stored  in  Eu3+:Y2SiOs  for  200  ps.  (c)  Image  reconstructed  under  the  same  condition  as 
that  of  (b)  but  with  weak  overall  intensity,  (d)  Image  reconstructed  5  minutes  after  the  completion  of 
recording.  In  (a)-(c),  the  data  pulse  length  was  1 1 .6  ps,  and  that  in  (d)  was  50  ps.  The  former  corresponds 
to  an  I/O  rate  of  300  megabits  per  second  while  the  latter  corresponds  to  67  megabits  per  second. 
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Figure  7.  (a)  Histogram  showing  the  intensity  distributions  of  0’s  and  I’s  that  were  retrieved  from  the 
time-domain  memory,  (b)  Normalized  histogram  obtained  after  applying  the  local  threshold  technique. 
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Poisson  distributions  to  accurately  determine  the  state  of  the  bit,  provided  that  the  intensity  within 
a  page  varies  smoothly.  We  adopted  the  following  formula  for  the  threshold 


lthr  =Io  +  P^SJo  +  4  (71  “ P^h)  ~ (h  +  P-f1 0 )] 


where  10  and  I ,  are  the  averages  of  0’s  and  l’s,  q  is  a  weighting  factor  chosen  empirically  to  be 
1/3,  and  p  is  a  measure  of  how  many  standard  deviations  above  I0  and  below  I,  we  can  place  the 
threshold.  Here,  we  use  p- 11.  The  optimum  numbers  used  in  averaging  were  found  to  be  four 
forO’s  and  eight  for  l’s. 

Obviously,  we  need  to  provide  initial  I0  and  /,  to  jump  start  the  analysis  of  each  page.  This 
is  accomplished  by  placing  4  known  0-bits  permanently  at  the  comers  of  a  4x4  region  in  the 
center  of  the  SLM  where  the  signal-to-noise  rate  is  largest.  We  average  the  4  bits  to  estimate  I0.  A 
crude  estimate  for  1,  can  be  obtained  by  averaging  the  intensities  of  the  remaining  twelve  unknown 
bits.  The  threshold  for  the  block  is  then  calculated  from  the  above  equation  by  assuming  the 
block  intensity  /Wdt=  (I0  +  /,)/ 2  which  implies  a  equal  number  of  l’s  and  0’s  in  the  block.  This 
approach  yields  a  BER  of  10‘6  for  the  12-bit  central  region,  or  equivalently  a  10'5  probability  of 
having  an  error  bit  per  page. 

We  then  continue  this  analysis  process  by  spiraling  outward  from  the  center  of  a  page 
towards  the  edges  with  the  intensity  decreasing  gradually  as  the  analysis  proceeds.  The 
neighboring  l’s  and  0’s  are  stored  in  two  separate  register  arrays  and  updated  in  a  first-in-first-out 
manner  to  obtain  their  running  averages.  Thus,  each  data  bit  is  determined  in  real  time  with  the  use 
of  the  equation  above.  In  comparison  to  the  existing  schemes  for  reading  page-formatted  digital 
data  [14],  this  automatic  threshold  calibration  technique  uses  no  coding  bits  and  requires  no 
mathematical  iteration.  These  advantages  lead  to  a  significant  increase  in  the  readout  rate,  ideal  for 
high-speed  applications. 

System  performance  and  projections.  Using  this  new  approach,  we  repeated  the 
above  experiment  and  obtained  a  histogram  for  the  same  100  pages  used  earlier  (Figure  7b).  This 
histogram  was  normalized  to  the  respective  threshold  of  each  data  bit  and  it  shows  two  clearly 
separate  distribution  curves.  The  separation  is  found  to  be  approximately  40%  of  the  full  width  of 
the  distribution  for  0’s.  This  large  separation  suggests  an  effective  increase  of  the  system’s 
dynamic  range  by  -40%.  The  raw  BER  is  estimated  to  be  ~10'7.  We  further  verified  this  number 
by  repeating  the  experiment  over  10  times  and  found  that  the  two  distribution  curves  in  the 
resulting  histogram  remain  separate. 

In  obtaining  the  above  results,  we  chose  the  data  pulse  to  be  1 1.2  ps  long,  which 
translates  to  an  I/O  rate  of  300  megabits  per  second.  The  100  pages  were  wavelength-division 
multiplexed  at  a  single  storage  cell  approximately  1.0  mm  x  1.0  mm  x  7.0  mm  at  a  rate  of  30 
frames  per  second  with  a  total  recording  time  of  -3.3  seconds.  All  the  pages  used  in  the 
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histogram  were  obtained  from  the  same  experiment.  For  convenience,  each  data  page  used  for 
BER  analysis  was  retrieved  after  being  transferred  to  the  memory  for  200  ps.  Using  short 
storage  time  prevents  undesired  page-to-page  fluctuation  caused  by  slow  frequency  drift  of  our 
laser.  The  echo  (or  “diffraction”)  efficiency  for  the  short  term  storage  was  measured  to  be 
-5  x  10'3. 

According  to  our  BER  measurement,  our  current  storage  system  is  capable  of  recording 
digital  data  close  to  10  megabits  without  generating  an  error  bit.  To  demonstrate  it,  we  performed 
an  experiment  in  which  we  stored  a  high  resolution  black-white  photograph  of  Albert  Einstein 
(Figure  8a).  This  digitized  photograph  contains  106  pages  totaling  356  kbits  of  digital  data. 
Figure  8b  shows  the  reconstructed  photograph  after  being  copied  to  and  read  back  from  the 
memory.  The  entire  photograph  was  stored  at  a  single  spatial  location  in  Eu3+:Y2Si05  using  WDM 
with  a  peak  recording  speed  of  300  megabits  per  second  and  no  error  correcting  code  of  any  kind 
was  used.  The  total  material  bandwidth  used  was  only  -85  MHz. 

The  present  work  represents  the  first  experimental  demonstration  of  parallel  digital  data 
storage  utilizing  the  time-domain  approach.  It  opens  up  the  possibility  of  developing  the 
technology  into  a  high-performance  page-oriented  dynamic  optical  memory.  As  the  technology 
for  parallel  data  access  becomes  more  mature,  this  time-domain  holographic  approach  along  with 
other  optical  holographic  storage  techniques  (e.g.,  the  photorefractive  memory)  may  be  able  to 
fulfill  their  long  awaited  promise  for  high-speed  high-density  data  storage. 

However,  a  number  of  issues  need  to  be  addressed  before  this  technology  can  be 
considered  commercially  viable.  One  of  them  is  the  low  operating  temperature  currently  required. 
We  hope  that  the  present  work  will  generate  sufficient  interests  in  the  scientific  community  so  that 
these  important  but  difficult  issues  can  be  addressed  and  hopefully  resolved  promptly. 
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Original 


Reconstructed 


(a) 


(b) 


Figure  8.  (a)  Digitized  photograph  displayed  using  its  original  binary  data,  (b)  Same  photograph  obtained 
after  the  binary  data  being  transferred  to  and  then  recalled  back  from  the  time-domain  memory. 
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CONCLUSIONS  AND  RECOMMENDATION 


In  conclusion,  we  have  demonstrated  the  practical  viability  of  the  stimulated  echo  approach 
using  rare-earth  doped  material  to  buffer  optical  data.  This  was  demonstrated  in  bit  serial  and  bit 
parallel  approach.  In  addition,  we  have  identified  new  materials  which  allow  the  storage  of  data  in 
the  1.5  itm  optical  range.  Optical  switching  of  data  packets  in  the  next  generation  optical  networks 
is  therefore  feasible. 
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APPENDIX  A 


OPTIMIZATION  OF  TIME-DOMAIN  STORAGE  DENSITY  IN  THE 
PRESENCE  OF  EXCITATION-INDUCED  SPECTRAL  DIFFUSION 


Optimization  of  time-domain  storage 
density  in  the  presence  of  excitation-induced 
spectral  diffusion 


X.  A.  Shen  and  R.  Kachru 


We  present  a  technique  to  reduce  the  loss  of  storage  capacity  in  time-domain  optical  memory  in  the 
presence  of  excitation- induced  spectral  diffusion.  We  applied  this  technique  to  a  Eu3":Y2Si05-based 
memory  that  utilizes  frequency  chirp  for  recording-retrieval  and  successfully  stored  860  bits  of  data  in 
a  40-MHz  data  channel.  In  comparison  with  the  conventional  chirp  method  under  identical  experimen¬ 
tal  conditions,  the  present  technique  offers  an  increase  of  a  factor  of  ~2  in  storage  capacity  with  no 
measurable  loss  of  signal  fidelity.  The  results  are  interpreted  in  terms  of  the  phase  error  accumulated 
during  dephasing-rephasing  and  spectral  diffusion  on  excited-state  relaxation.  ©  1997  Optical  Society 
of  America 

Key  words:  Time- domain  optical  memory,  photon  echoes,  frequency  chirp,  excitation-induced  spec¬ 
tral  diffusion. 


Since  proposed  in  the  early  1980’s,1  coherent  time- 
domain  optical  memory  (TDOM)  has  emerged  as  a 
strong  candidate  for  next-generation  storage  devices. 
In  TDOM,  data  can  be  stored  and  recalled  rapidly  in 
either  bit-serial  or  bit-parallel  formats,  providing  un¬ 
precedented  flexibility  for  diversified  applications 
varying  from  optical  packet  switching  to  high-speed 
image  processing.  Single-channel  storage-retrieval 
of  serial  data  at  40  Mbits/s,  parallel  recording  at  an 
input- output  speed  of  —20  (xs/frame,  as  well  as  a 
storage  density  of  ~8  Gbits/in.2  (~~8  Gbits/6.45  cm2) 
have  been  demonstrated  with  cw  lasers  of  a  modest 
peak  power.2-4 

While  the  theoretical  limits  of  the  bandwidth  and 
storage  capacity  are  substantially  high  for  TDOM  as 
a  result  of  the  presence  of  the  frequency  dimension 
for  recording  (typically  5-10  GHz  and  106-10 '  bits/ 
focal  volume  for  rare-earth-doped  solids),  reaching 
these  limits  remains  a  challenge.  One  of  the  obsta¬ 
cles  is  the  presence  of  excitation-induced  spectral  dif¬ 
fusion  (EISD)5-9  that  results  from  interactions  among 
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active  ions  during  information  recording-retrieval. 
This  process  leads  to  a  broadening  of  the  recorded 
spectral  gratings,  which  in  turn  reduces  the  useable 
storage  capacity  of  the  memory. 

Here  we  present  a  technique  to  maximize  the  stor¬ 
age  capacity  of  TDOM  in  the  presence  of  EISD.  We 
choose  frequency-chirp  recording  as  an  example  and 
show  that,  with  the  proposed  technique,  an  increase 
of  a  factor  of  —2  in  the  storage  density  can  be 
achieved  in  a  Eu3+:Y2Si05-based  memory. 

In  TDOM,  information  is  stored  in  a  recording  me¬ 
dium  by  a  write  pulse  Ew ,  which  often  precedes  the 
data  pulse  by  a  time  t  and  serves  as  a  temporal 
reference  for  data  retrieval.  To  retrieve  the  infor¬ 
mation  T  (storage  time)  after  the  arrival  of  the  data 
pulse,  the  recording  medium  is  reexposed  to  a  read 
pulse  Er,  which  induces  coherent  radiation  Ee 
known  as  stimulated-photon  echo.  Ee  occurs  also  at 
a  time  t  after  the  arrival  of  ER  and  is  related  to  the 
three  inputs  as  follows: 

Ee{v)  a  Ew*(v)ER(v)ED(v)f  (1) 

where  Et  is  the  Fourier  transform  of  Et. 

In  TDOM,  EISD  manifests  itself  as  a  loss  of  echo 
intensity,  but  in  two  distinct  manners6-9:  a  rapid 
decrease  of  instantaneous  echo  power  as  t  increases, 
and  a  decrease  of  overall  echo  intensity  as  T  in¬ 
creases.  The  former  is  known  as  excitation-induced 
excess  dephasing,  which  is  found  to  result  from  phase 
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Fig.  1.  (a)  Power  spectrum  of  a  linearly  chirped  pulse  with  a  rectangular  envelope  {solid  curve)  and  (b)  the  power  spectrum  of  an 

identically  chirped  pulse  (solid  curve)  with  a  truncated  Gaussian  envelope  (the  insert).  Here,  A  is  the  chirp  width  and  Tp±  =  256.  The 
dashed  curves  in  both  (a)  and  (b)  are  the  same  spectrum  of  a  brief  rectangular  pulse  of  length  t  =  2/A. 


errors  accumulated  by  active  ions  during  data  record’ 
ing  (0  — »  t,  dephasing  time)  and  retrieval  (T  +  t  — ► 
T  +  2t,  rephasing  time).6'10  These  errors  are  gen¬ 
erated  by  interaction  with  neighboring  ions  undergo¬ 
ing  a  change  of  electronic  states  on  laser  excitation. 
These  ions  modify  the  electric  field  seen  by  an  active 
ion,  giving  rise  to  a  shift  of  its  transition  frequency. 
The  larger  the  time  t,  the  higher  the  probability  for 
this  ion  to  accumulate  a  phase  error. 

When  t  — >  0,  excess  dephasing  vanishes,  and  spec¬ 
tral  diffusion  becomes  dominant.  Depending  on  the 
choice  of  T,  a  number  of  diffusion  mechanisms  may  be 
present.11 12  For  values  of  T  comparable  with  the 
excited-state  lifetime  Tx,  diffusion  from  excited-state 
ion  relaxation  often  dominates.  This  process  results 
in  a  broadening  of  the  spectral  grating  generated  by 
Ew  and  ED ,  which  leads  to  a  loss  of  echo  intensity, 
hence  a  loss  of  storage  capacity. 

EISD  is  determined  by  the  concentration  of 
excited-state  ions.  For  a  spectrally  and  spatially 
uniform  ion  distribution,  EISD  is  then  determined  by 
the  energy  of  the  excitation  pulses.  The  lower  the 
pulse  energy,  the  weaker  the  EISD.  Thus,  prevent¬ 
ing  EISD-mediated  memory-capacity  loss  reduces  to 
finding  laser  pulses  that  offer  a  maximum  Ee  at  a 
fixed  pulse  energy. 

For  convenience,  we  assume  the  following:  (a) 
The  data  stream  to  be  stored  consists  of  only  one 
short  rectangular  pulse  ED  of  a  length  t  and  an  am¬ 
plitude  Ed o,  and  (b)  Ew  —  ER.  To  make  our  analysis 
independent  of  EISD,  we  compare  different  recording 
methods  under  identical  energies  for  Ew  and  ER. 

We  first  examine  the  conventional  chirp  meth¬ 
od13’14  for  which  we  obtain  Ew  and  ER  by  linearly 
sweeping  the  carrier  frequency  of  a  fixed-power  laser 
pulse.  Figure  1(a)  shows  the  power  spectrum  (the 
solid  curve)  of  such  a  pulse  with  A Tp  —  256,  where  A 
is  the  chirp  width  and  Tp  is  the  pulse  length.  The 
dashed  curve  in  Fig.  1(a)  is  the  power  spectrum  of  the 
data  pulse  defined  above.  We  chose  A  —  2/t  (the  full 
width  of  the  data  pulse  at  the  first  minima).  For 
minimizing  the  EISD,  A  must  be  equal  to  —2/t  be¬ 
cause  at  a  given  chirp  rate  a  larger  A  means  more 
energy  is  deposited  outside  of  the  data  bandwidth. 


and  those  energies  contribute  only  weakly  to  Ee  but 
equally  to  EISD. 

Even  under  the  condition  of  A  =  2/t,  the  conven¬ 
tional  chirp  method  deposits  an  equal  amount  of  en¬ 
ergy  at  the  band  edges,  where  the  data  signal  nearly 
vanishes  [see  Fig.  1(a)].  This  implies  that  the  con¬ 
ventional  method  does  not  use  the  energy  of  Ew  and 
Er  efficiently.  If  EISD  were  present,  it  would  cer¬ 
tainly  result  in  a  small  echo  signal  and  consequently 
a  low  storage  density. 

One  approach  that  utilizes  the  pulse  energy  effi¬ 
ciently  is  to  chirp  Ew  and  ER  in  such  a  way  that 
their  power  spectra  resemble  closely  that  of  the 
data  pulse.  Since  the  desired  spectra  are  identical 
to  that  of  the  data  pulse,  we  can  estimate  Ee  gen¬ 
erated  by  such  chirped  pulses  by  using  simple,  brief 
Ew  and  ER  pulses  that  are  temporally  identical  to 
Ed.  For  a  brief  Ew  {=ER)  of  a  length  t  and  ampli¬ 
tude  E0,  the  resulting  peak-echo  amplitude  EE0brief  * 
3E02ED0t2/4  with  a  full  width  at  half-maximum  of 
-1.26t. 

We  can  compare  the  above  amplitude  to  that  ob¬ 
tained  with  the  conventional  chirp  method.  Notice 
that,  for  a  rectangular  chirped  Ew  (here  A  =  2/t)  to 
have  the  same  energy  as  the  above  brief  pulses,  its 
spectrum  peak  \Ew\max2  must  satisfy  \Ew\max2h  » 
E2 t.  From  relation  (1)  we  obtain  a  peak-echo  am¬ 
plitude  Ee 0rect  «  E02ED0t2/2,  with  a  width  of  — t, 
where  the  superscript  rect  denotes  rectangular. 
Therefore,  we  have  a  ratio  of  EEOhrie{/EEOrect  =  3/2 
and  a  corresponding  intensity  ratio  of  —9/4. 

The  above  analysis  states  that,  in  the  presence  of 
EISD,  an  increase  of  a  factor  of  —2.25  in  the  echo 
power  can  be  achieved  by  redistribution  of  the  energy 
of  Ew  and  ER  in  such  a  way  that  their  power  spectra 
are  identical  to  that  of  the  data  pulse  to  be  stored. 
In  TDORTs  utilizing  frequency  chirp  for  recording- 
retrieval,  this  can  be  accomplished  by  appropriate 
modulation  of  the  amplitudes  of  Ew  and  ER  so  that 
the  energy  at  the  band  edges  (which  corresponds  to 
the  early  and  later  parts  of  Ew  and  ER)  is  reduced. 

Figure  1(b)  shows  the  power  spectrum  obtained  by 
replacement  of  the  conventional  chirped  pulse  [Fig. 
1(a)]  with  a  chirped  pulse  whose  envelope  resembles 
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Fig.  2.  Autocorrelations  of  the  same  pulses  shown  in  Figs.  1(a) 
(dashed  curve)  and  Kb)  (solid  curve).  Here,  A  is  the  chirp  width. 

a  truncated  Gaussian  curve  of  the  same  length  [Fig. 
1(b)].  Here,  the  Gaussian  width  (full  width  at  1/e) 
was  chosen  to  be  — 0.78Tp.  Calculations  show  that 
the  echo  intensity  obtained  with  such  pulses  is  higher 
by  a  factor  of  —2.0  than  that  obtained  with  conven¬ 
tional  chirped  pulses  of  the  same  energy  and  that  the 
width  of  the  retrieved  pulse  is  —  1.26t,  which  is  in 
good  agreement  with  the  estimates  obtained  with 
brief  Ew  and  ER  pulses. 

In  digital  data  storage,  the  broadening  of  a  sub¬ 
pulse  in  a  data  stream  leads  to  an  increase  of  the 
signal  level  between  two  adjacent  pulses.  If  the 
peak  of  a  subpulse  represents  1  and  the  absence  of  a 
pulse  represents  0,  this  broadening  translates  to  a 
broadening  of  the  intensity  distribution  of  0’s.  Care¬ 
ful  examination  of  F~l{Ew*(v)ER(v)}  (the  inverse 
Fourier  transform  representing  the  correlation  of  ER 
with  Ew)  for  the  above  two  methods  shows  that  the 
conventional  chirp,  although  it  offers  a  sharper  cor¬ 
relation  peak,  generates  large  sidelobes  (the  dashed 
curve  in  Fig.  2),  which  cause  a  fluctuation  of  the  echo 
amplitude  and  consequently  a  broadening  of  the  dis¬ 
tribution  of  l’s.  The  proposed  chirp  method,  on  the 
other  hand,  has  a  broadened  correlation  peak  but 
negligible  sidelobes  (the  solid  curve  in  Fig.  2),  the 
latter  reducing  the  intensity  fluctuation  of  l’s.  Cal¬ 
culations  show  that  the  sum  of  the  two  distribution 
widths  (a  measure  of  the  system’s  bit-error  rate)  is 
actually  reduced  by  —27%  with  the  use  of  the  pro¬ 
posed  method. 

In  the  implementation  of  TDOM,  Ew  and  ER  are 
often  obtained  by  modulation  of  a  cw  laser 
externally.2'4  The  maximum  achievable  peak 
power  is  thus  fixed.  We  have  shown  that,  if  EISD 
were  present,  using  the  conventional  chirp  method 
with  a  reduced  peak  power  would  result  in  a  loss  of  a 
factor  of  —2  in  the  signal  intensity,  which  in  turn 
leads  to  a  loss  of  storage  capacity.  The  only  question 
that  remains  to  be  answered  is  which  of  the  two  chirp 
methods  would  offer  a  higher  storage  density  if  both 
were  operated  at  the  same  peak  power.  Obviously, 
the  conventional  method  would  give  a  larger  echo 
signal  if  EISD  were  absent  because  more  energy 
would  be  used.  If  EISD  were  present,  the  answer 
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Fig.  3.  Data  retrieved  from  a  Eu"3:Y2Si05-based  memory  by  use 
of  three  different  pairs  of  Ew  and  ER :  (a)  linearly  chirped  Ew  and 
Er  with  a  rectangular  envelope,  (b)  identically  chirped  Ew  and  ER 
with  a  truncated  Gaussian  envelope  (full  width  at  1/e  =  0.7 8Tp), 
and  (c)  identically  chirped  Ew  and  ER  with  an  envelope  correspond¬ 
ing  to  the  first  half  of  a  sine  function.  In  all  three  cases,  the  peak 
laser  power,  Tpt  and  A  remain  the  same,  and  they  are  -100  mW, 
6.4  p. s,  and  40  MHz,  respectively. 


would  depend  on  the  storage  material  and  would 
have  to  be  determined  experimentally. 

We  performed  an  experiment  on  the  7F0-5D0  tran¬ 
sition  in  Eu3+:Y2Si05  to  address  this  question.  We 
chose  Eu3+:Y2Si05  because  in  this  material  diffu¬ 
sions  other  than  EISD  were  found  to  be  insignificant 
below  4  K,8  and  any  loss  of  echo  intensity  could  be 
attributed  primarily  to  EISD.  The  experimental  ap¬ 
paratus  and  conditions  were  identical  to  those  de¬ 
scribed  in  Ref.  15. 

We  stored  identical  data  streams  that  contain  860 
bits  of  data  using  three  different  chirp  methods:  (a) 
the  conventional  chirp  method,  (b)  chirped  Ew  and  ER 
with  the  truncated  Gaussian  amplitude  envelope  dis¬ 
cussed  above,  and  (c)  chirped  Ew  and  ER  with  an 
envelope  equal  to  the  first  half  of  the  sine  function. 
In  all  three  cases,  the  peak  power  of  Ew  and  ER ,  Tp, 
A,  the  laser  beam  waist,  the  separation  between  Ffw 
and  ED)  and  the  separation  between  ED  and  ER  are 
chosen  to  be  identical,  and  they  are,  respectively, 
—  100  mW,  6.4  (jls,  40  MHz,  —100  pm  (FWFM),  18  jjls, 
and  240  |xs.  The  length  of  each  data  subpulse  was 
-50  ns,  giving  rise  to  a  data  length  of  —90  \xs  and  a 
bandwidth  equal  to  A.  The  pulse  energy  used  in 
method  (a)  was  —2.1  times  larger  than  that  of  meth¬ 
ods  (b)  and  (c). 

Figures  3(a)-3(c)  show  the  experimental  results  ob¬ 
tained  with  the  above  three  methods,  respectively. 
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The  peak-echo  powers  in  all  three  cases  are  nearly 
identical.  The  conventional  chirp  method  (a)  yields 
a  faster  decay  in  comparison  with  the  other  methods, 
which  show  decay  rates  comparable  with  each  other. 
If  the  signals  at  the  ends  of  the  data  strings  shown  in 
Figs.  3(b)  and  3(c)  correspond  to  the  minimum  resolv¬ 
able  signal  level  of  our  system,  then  the  fast  signal 
decay  shown  in  Fig.  3(a)  translates  to  a  loss  of  a  factor 
of  —2  in  data  capacity.  The  exact  gain  of  the  storage 
capacity  requires  careful  examination  of  the  fidelity 
of  each  data  bit.  However,  visual  inspection  of  the 
echo  data  shows  no  apparent  degradation  of  signal 
fidelity  with  the  use  of  the  proposed  technique  (see 
the  inserts  in  Figs.  3). 

Plotting  the  results  on  a  logarithmic  scale  shows  that 
the  decays  closely  follow  the  function  exp(-rex^)1  with 
~~  1/(33  fis)  for  Fig.  3(a)  and  T ^  -  1/(65  jxs)  for  Figs.  3(b) 
and  3(c).  Extrapolation  of  the  data  to  the  center  of  the 
read  pulse  (i.e.,  t  0)  yields  an  echo-intensity  ratio  be¬ 
tween  the  results  of  Figs.  3(a)  and  3(b)  [or  Figs.  3(a)  and 
3(c)]  of  —1.2,  much  smaller  than  the  value  of  2.2  calcu¬ 
lated  in  the  absence  of  EISD.  This  difference  can  be 
attributed  to  partial  relaxation  of  the  excited-state  ions 
generated  by  the  additional  (-51%)  pulse  energy  used  in 
the  conventional  chirp  method.  The  ratio  between  the 
two  decay  constants  [i.e.,  (65  jjls)/(33  (jls)  =  1.97]  is  a 
direct  measure  of  the  gain  in  storage  capacity  obtained 
with  the  proposed  technique. 

In  conclusion,  we  have  presented  a  novel  technique 
to  increase  the  useable  storage  capacity  of  a  TDOM  in 
the  presence  of  EISD.  This  technique  has  been  tested 
experimentally  in  Eu3+:Y2Si05  and  an  increase  of  a 
factor  of  —2  in  storage  capacity  has  been  achieved. 

We  express  our  appreciation  to  John  W.  Perry  for 
the  numerical  simulation.  This  research  was  sup¬ 
ported  by  the  Defense  Advanced  Research  Projects 
Agency  under  contract  F-49620-95-C-0077. 
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We  report  what  is  to  our  knowledge  the  first  experimental  observation  of  two-  and  three-pulse  echoes  in 
Er3":YAG  at  1.527  ^m.  The  experiment  was  performed  on  the  transition  from  the  AI\sn  ground  state  to  the 
lowest  Stark  level  of  the  4/13, 2  excited  state.  The  corresponding  homogeneous  dephasing  time  was  measured  to 
be  approximately  50  ns  at  1  and  7  K.  In  addition,  three-pulse  echo  decay  indicated  a  lifetime  of  the  spectral 
grating  of  —30  fi s.  Implications  of  these  results  for  optical  data  storage  and  optical  signal  processing  are 
discussed.  ©  1997  Optical  Society  of  America 


Photon  echo  in  solids  provides  not  only  important  spec¬ 
troscopic  information  about  the  dynamic  interactions 
among  active  ions  and  their  host  materials  but  also 
a  powerful  means  of  storing  and  processing  temporal 
signals  of  various  kinds  at  optical  wavelengths;  the 
latter  feature  led  to  the  proposal  and  successful  ex¬ 
perimental  demonstration  of  coherent  time-domain 
optical  memory  (TDOM).1-5  Unlike  other  optical 
memories,  TDOM  is  a  high-speed  signal  processor  by 
nature.  Temporal  data  strings  stored  at  different 
times  can  be  correlated  or  convoluted  simultaneously 
with  an  incoming  data  string  in  the  memory  at  band- 
widths  of  up  to  several  gigahertz,  with  no  additional 
processing  time  required.6-8  This  rapid  in-memory 
processing  feature  combined  with  the  demonstrated 
high  storage  capacity  also  makes  TDOM  a  superior 
candidate  for  applications  such  as  high-speed  optical 
packet  switching,  context-based  retrieval,  matched 
filtering,  data  encryption,  and  pattern  recognition. 

In  spite  of  the  rapid  progress  in  the  development 
of  optical  memories  and  all-optical  processors,  to  our 
knowledge  there  exists  no  model  directly  operating  at 
optical  communication  wavelengths.  The  current  op¬ 
tical  devices  require  electrical -optical  signal  conver¬ 
sions,  which  certainly  reduce  the  data-communication 
speed.  Meeting  the  rapidly  increasing  demand  for 
high-speed  and  high-capacity  optical  communications 
has  made  it  necessary  to  eliminate  such  electrical - 
optical  conversions.  However,  it  remains  a  challenge 
to  find  a  suitable  medium  for  storing  and  processing 
optical  signals  at  the  infrared  wavelengths  of  fiber 
communication. 

In  this  Letter  we  explore  the  feasibility  of  developing 
a  TDOM  that  can  be  operated  directly  at  the  communi¬ 
cation  wavelengths.  We  choose  the  47i5/2-47i3/2  transi¬ 
tion  in  Er3+:YAG  as  a  candidate  for  the  present  study 
and  examine  some  of  its  spectroscopic  properties  that 
are  relevant  to  time-domain  applications. 

The  4f  15/2  13/2  transition  of  triply  ionized  Er  is 
widely  used  in  lasers  and  optical  amplifiers.  This 
transition,  occurring  at  wavelengths  near  1.5  yum,  suit¬ 
able  for  fiber  communications,  has  been  studied  in 
glass  fiber9-11  and  crystal.12  These  studies  focused 
on  either  the  absorption  and  the  fluorescence  of  the 
transition10-12  or  on  picosecond  photon  echo  in  Er- 
doped  fiber.9  Little  was  known  about  the  homoge¬ 


neous  dephasing  and  spectral  diffusion  in  Er3* -doped 
crystal;  these  properties  are  critical  to  the  application 
of  TDOM. 

We  first  performed  an  absorption  measurement  of 
Er3*:YAG  to  resolve  the  4/i5/2-47i3/2  transition  by 
using  a  tunable  pulsed  dye  laser  (Quanta  Ray  PDL-2) 
that  was  Raman  shifted  in  a  H2  gas  cell.  The  laser 
was  operated  at  a  repetition  rate  of  10  Hz  with  a  pulse 
width  of  —7  ns  and  pulse  energy  of  —2  mJ.  The 
accuracy  of  the  laser  wavelength  was  estimated  to  be 
0.1  nm.  The  linewidth  of  the  laser  was  estimated  to 
be  -4  GHz.  The  sample  of  Er3*:YAG  has  an  Er3* 
concentration  of  —0.1%  and  was  placed  either  in  a 
closed-cycle  refrigerator  at  —7  K  or  in  a  liquid-helium 
Dewar  at— IK. 

The  Er  ion  has  an  odd  number  of  f  electrons, 
giving  rise  to  Kramers  degeneracy,  which  can  be  lifted 
only  by  an  external  magnetic  field.  We  resolved  six 
lines  associated  with  transitions  from  the  lowest  Stark 
level  in  the  4/is/2  state  to  six  Stark  levels  of  the 
4fi3/2  state  (Fig.  1).  Their  wave  numbers  from  right 
to  left  in  Fig.  1  are  6549,  6599,  6606,  6786,  6806, 
and  6886  cm-1.  The  absorption  line  at  6527  cm'1 
at  —7  K  is  a  hot  band  from  the  second-lowest  Stark 
level  (—22  cm-1)  in  the  4/is/2  state  to  the  lowest 
level  of  the  state.  At  higher  temperature,  we 

also  observed  a  hot  band  from  higher  stark  levels 
of  the  4/i3/2  state.  The  energy-level  diagram  of  4/is/2 
and  47i3/2  that  corresponds  to  our  data  is  shown  in 
Fig.  2.  Our  results  are  in  general  consistent  with 
research  on  Er3*-YAG  by  Ashurov  et  a/.,13  who  also 
observed  only  six  lines  for  the  4/is/2-47i3/2  transition  at 
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Fig.  2.  Energy-level  diagram  of  Er3+:YAG  for  the  4/is/2 
and  47i3/2  state. 


4  K.  The  data  reported  by  Koningstein  and  Geusic14 
indicate  an  additional  absorption  line  —91  cm'1  above 
the  highest  line  in  our  data,  which  may  be  too  weak  for 
us  to  detect. 

The  rest  of  our  experiment  focuses  on  the  tran¬ 
sition  from  the  ground  state  of  4/i5/2  to  the  low¬ 
est  Stark  level  of  the  AI\%i2  state  (i.e.,  6549  cm"1,  or 
1.527  /xm).  First  the  Al\zn  excited-state  lifetime  was 
measured  with  a  pump -probe  method.  We  resonantly 
excited  Er  ions  to  that  level  with  an  infrared  dye 
laser.  We  then  probed  the  population  of  that  level 
by  applying  another  pulsed  dye  laser  to  excite  the 
Er  ions  from  the  AI\zn  state  to  the  4F7/2  state  (at 
20  516  cm"1)  and  detecting  the  fluorescence  from  the 
4F7/2-  By  delaying  the  probe  pulse  with  respect  to 
the  pump  pulse,  we  measured  the  47i3/2  state  lifetime 
to  be  —10  ms,  consistent  with  the  result  reported  by 
Li  et  aZ.15 

The  inhomogeneous  broadening  of  the  transition 
was  estimated  to  be  —20  GHz  from  the  measured 
transmission  spectrum  (Fig.  1).  We  measured  the  ho¬ 
mogeneous  linewidth  of  this  transition  by  recording 


two-pulse  echo  intensity  as  a  function  of  excitation 
pulse  separation  [shown  in  Fig.  3(a)] .  A  fit  of  the  data 
to  a  simple  exponential  function  gives  a  homogeneous 
dephasing  time  T2  of  —50  ns.  This  dephasing  time 
was  found  to  vary  little  between  1  and  7  K  under  a 
magnetic  field  of  less  than  100  G. 

In  our  two-pulse  echo  experiment  we  use  two  dye 
lasers  that  are  similar  to  the  one  described  for  the  ab¬ 
sorption  measurement.  The  peak  intensity  of  the  laser 
inside  the  crystal  is  —1  MW /cm2.  We  use  the  fluo¬ 
rescence  lifetime  of  the  47i3/2  excited  state  to  calculate 
the  laser  pulse  area  and  take  into  account  that  the 
laser  pulse  is  non-Fourier-transform  limited.16  The 
estimated  pulse  area  is  between  77/8  and  tt/4.  This 
area  of  <1  was  also  supported  by  the  experimental  ob¬ 
servation  that  the  echo  intensity  was  linear  to  the  first 
laser  pulse’s  intensity  and  quadratic  to  the  intensity 
of  the  second.  We  also  found  that  reducing  the  laser 
power  by  a  factor  of  2  had  no  effect  on  the  decay  rate  of 
the  echo. 

The  data  plotted  in  Fig.  3(a)  were  taken  with  an 
external  magnetic  field  of  —100  G  applied  along  the 
c  axis  of  the  crystal.  When  this  external  field  was 
increased  from  0  to  10  G,  the  echo  intensity  quickly 
increased  by  a  factor  of  3.  But  this  increase  became 
saturated  when  the  field  was  greater  than  30  G.  On 
the  other  hand,  the  dephasing  time  did  not  change 
when  we  increased  the  magnetic  field  to  100  G.  When 
a  magnetic  field  of  —3000  G  was  applied  to  the  crystal 
(perpendicular  to  its  c  axis),  we  observed  apparently 
nonexponential  decay  of  echo,  but  the  dephasing  time 
was  not  lengthened  either. 

One  of  the  important  parameters  governing  the 
performance  of  the  TDOM  is  the  highest  temperature 
at  which  the  storage  sample  can  be  operated.  To 
determine  this  temperature,  we  measured  two-pulse 
echo  intensity  at  different  sample  temperatures.  At 
a  fixed  excitation  pulse  separation  of  —35  ns,  the  echo 
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Fig.  3.  (a)  Two-pulse  echo  intensity  in  Er3":YAG  at  —7  K 
under  an  external  magnetic  field  of  100  G  as  a  function 
of  change  of  excitation  pulse  separation.  Open  circles 
are  data,  and  the  solid  line  is  a  simple  exponential  fit. 
(b)  Stimulated  photon-echo  intensity  in  Er3^:YAG  as  a 
function  of  the  separation  between  the  second  and  the 
third  excitation  pulses.  The  results  were  obtained  at  —1  K 
under  an  external  magnetic  field  of  3000  G. 
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intensity  gradually  decreased  by  only  a  factor  of  2  from 
7  to  1 1  K.  This  temperature  dependence  of  Er3+ :YAG 
is  considerably  weaker  than  that  of  other  rare-earth- 
doped  crystals,  except  for  Eu3+:Y2SiOs  at  its  'Fo-°Do 
transition.1' 

In  addition  to  performing  the  two-pulse  echo  experi¬ 
ments,  we  also  examined  three-pulse  (or  stimulated) 
photon  echo  in  Er3":YAG,  since  it  is  closely  related 
to  the  physical  processes  involved  in  TDOM  and  re¬ 
lated  optical  processors.  By  doing  so,  we  were  able 
to  estimate  the  potential  performance  of  Er3+:YAG  in 
TDOM  at  1.527  (i m,  for  example,  the  maximum  stor¬ 
age  time,  which  is  often  limited  by  spectral  diffusion 
resulting  from  ion -host  and  ion- ion  interactions.  We 
carried  out  the  three-pulse  stimulated  echo  experiment 
under  a  condition  very  similar  to  that  of  the  two-pulse 
echo  experiment.  We  fixed  the  separation  between  the 
first  two  excitation  pulses  at  20  ns  and  varied  that  be¬ 
tween  the  second  and  the  third  pulses.  The  results 
[shown  in  Fig.  3(b)]  show  a  complicated  decay  feature, 
with  an  initial  fast  decay  followed  by  a  slower  decay 
and  then  another  fast  decay.  At  a  delay  of  50  /ns  the 
echo  intensity  decreased  by  a  factor  of  10  from  that 
at  1  ps.  This  decay  is  much  faster  than  the  fluores¬ 
cence  lifetime  measured  in  this  study,  suggesting  the 
presence  of  rapid  spectral  diffusion  on  a  microsecond 
time  scale.  One  possible  cause  of  this  fast  diffusion  is 
the  strong  Er-Er  dipolar  interactions,  as  suggested  by 
Macfarlane  and  Shelby.18 

In  conclusion,  we  have  observed  for  what  is  believed 
to  be  the  first  time  two-  and  three-pulse  echoes  at 
1.527  pm  in  Er3^:YAG.  Although  the  fast  dephasing 
of  Er3+:YAG,  together  with  its  insensitivity  to  an  ex¬ 
ternal  magnetic  field  as  observed  in  our  experiments, 
is  less  desirable,  it  poses  no  significant  challenge  to 
the  implementation  of  Er3+:YAG  as  a  time-domain 
storage  medium  in  most  applications.  For  example, 
the  recently  developed  swept-carrier  technique  per¬ 
mits  the  storage  of  data  strings  that  are  much  longer 
than  the  dephasing  time  of  a  storage  material.  The 
lifetime  of  the  spectral  grating  measured  by  three- 
pulse  echoes,  although  it  is  short,  is  still  adequate  for 
applications  such  as  dynamic  optical  memory,  optical 
packet  buffering,  and  optical  signal  processing.  The 
results  reported  in  this  Letter  open  up  the  possibility 
of  developing  TDOM  and  TDOM-based  optical  proces¬ 
sors  operated  directly  at  communication  wavelengths. 
Further  improvement  of  material  properties,  such  as 
longer  dephasing  time  and  less  spectral  diffusion,  could 


be  accomplished,  e.g.,  by  choice  of  different  host  ma¬ 
terials,  reducing  Er  concentrations,  and  (or)  applying 
larger  external  magnetic  fields.  These  efforts  have 
shown  success  in  other  rare-earth  ions  such  as  Eu  and 
Pr.  However,  they  may  not  be  applicable  to  Er  ions. 

This  research  was  supported  by  the  Defense  Ad¬ 
vanced  Research  Projects  Agency  under  contract 
F-49620-95-C-0077.  Useful  discussions  with  Robert 
Leheny  and  Alan  Craig  are  greatly  appreciated. 
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APPENDIX  C 


TIME-DOMAIN  HOLOGRAPHIC  DIGITAL  MEMORY 


Time-Domain  Holographic  Digital  Memory 

Xiao  A.  Shen,  An-Dien  Nguyen,  John  W.  Perry,* 

David  L.  Huestis,  Ravinder  Kachrut 

An  optical  storage  technique  based  on  time-domain  holography  for  the  rapid  recording 
and  readout  of  page-formatted  digital  data  is  demonstrated.  Storage  of  356  kilobits  of 
data  was  achieved  at  a  single  spatial  location  in  a  rare-earth-doped  crystal.  The  digital 
data  were  recorded  and  accurately  retrieved  at  a  peak  rate  of  300  megabits  per  second 
without  the  use  of  error-correcting  codes.  The  system’s  raw  bit  error  rate  is  about  1 0“7. 
This  low  bit  error  rate  was  achieved  by  a  detection  scheme  for  extraction  of  binary  data. 
These  results  have  implications  for  dynamic  optical  memory. 


Holographic  memories,  once  considered  a 
storage  device  of  the  distant  future,  have 
come  much  closer  to  reality.  Instruments 
that  make  use  of  holographic  approaches  to 
data  storage  have  been  developed  and  tested 
for  applications  such  as  the  identification  of 
fingerprints  ( 1 ).  Use  of  spatial  holography  to 
store  page -formatted  digital  optical  data  in 
computer-based  applications  has  also  been 
explored,  and  the  results  are  promising  (2- 
4).  The  memories  are  believed  to  be  superior 
to  existing  technologies,  with  key  features 
including  ultrahigh  storage  density,  rapid 
data  transfer,  short  data  access  times,  and 
exceptional  reliability. 

Although  much  progress  has  been 
achieved  in  some  key  areas  of  holographic 
data  storage,  several  major  technical  imped¬ 
iments  remain,  restricting  its  application  in  a 
broad  area,  as  was  promised  initially.  The 
most  noticeable  obstacle  has  been  the  slow 
recording  rate  of  conventional  holographic 
techniques,  which  record  spatial  interference 
patterns  generated  by  two  temporally  co¬ 
existing  laser  beams.  The  recording  requires, 
in  addition  to  laser  excitation,  subsequent 
electron  diffusion  in  storage  materials,  as  in 
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the  case  of  photorefractive  memory,  for  ex¬ 
ample.  Such  processes  limit  photorefractive 
technology  to  applications  that  are  not  input 
intensive,  such  as  memory  that  is  written 
only  once.  Success  in  developing  high-per¬ 
formance  input-output  (I/0)-intensive  ho¬ 
lographic  memories,  such  as  page-oriented 
dynamic  optical  random  access  memory,  will 
require  the  development  of  much  higher  re¬ 
cording  speeds  while  maintaining  high  stor¬ 
age  densities  and  low  bit  error  rates,  the 
latter  being  a  benchmark  for  the  quality  of 
stored  data. 

Time-domain  holographic  memory  has 
emerged  as  an  excellent  candidate  for  such 
I/O  intensive  applications  (5,  6).  This  mem¬ 
ory  not  only  has  features  common  to  con¬ 
ventional  spatial  holographic  recording, 
such  as  high  capacity  and  high  degree  of 
parallelism,  but  also  some  that  are  comple¬ 
mentary  to  conventional  holography,  in¬ 
cluding  fast  page  recording  speed  and  modest 
laser  power  requirement.  Recent  experi¬ 
ments  (6)  show  that  an  exposure  time  of 
several  microseconds  and  a  peak  laser  power 
of  less  than  200  mW  are  sufficient  to  record 
a  high-resolution  binary  image  with  the 
time-domain  technique. 

Historically,  the  time-domain  holo¬ 
graphic  approach  evolved  from  an  optical 
transient  phenomenon  known  as  stimulated 
photon  echoes  (7).  It  is  a  time-domain  or, 
equivalently,  a  spectral -domain  version  of 
holographic  recording.  Like  conventional 
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Fig.  1.  Time-domain  holo¬ 
graphic  recording  scheme. 
The  insert  shows  the  order 
in  which  the  write,  data, 
read,  and  reconstructed 
pulses  occur. 


Read  pulse 


holography,  chis  technique  stores  data  as 
holograms  by  using  a  reference  beam  and  a 
data  beam.  These  two  beams  interfere  at  a 
recording  medium,  causing  a  change  of  the 
optical  properties  of  the  medium.  Unlike  in 
spatial  holography,  the  reference  and  data 
beams  appear  in  the  form  of  short  laser 
pukes  and  do  not  overlap  temporally.  In 
most  recording  schemes  developed  to  date, 
the  reference  pulse  precedes  the  data  pulse 
by  several  nanoseconds  to  several  microsec¬ 
onds  to  provide,  in  addition  to  a  spatial 
reference,  a  temporal  reference  for  the  re¬ 
corded  holograms  (Fig.  1).  The  interference 
generated  by  these  two  temporally  distinct 
laser  pulses  occurs  in  both  spectral  and  spa¬ 
tial  domains,  giving  rise  to  four-dimensional 
holograms  (8). 

The  spectral  interference  of  the  two 
pukes  is  recorded  in  a  storage  medium 
[which  is  often  a  rare -earth- doped  solid 
such  as  Eu3TY2SiC>5  or  PHTYAG  (YAG, 
yttrium-aluminum-gamet)]  by  laser  excita¬ 
tion.  Ions  such  as  Eu3+  and  Pr3*  that  are 
responsible  for  data  recording  are  resonant¬ 
ly  excited  by  the  laser  pukes  to  form  popu¬ 
lation  gratings  that  resemble  the  Fourier 
spectra  of  the  pulse  pair  (Fig.  2).  These 
gratings  in  a  carefully  chosen  material  can 
persist  for  several  hours  to  perhaps  several 
days  at  cryogenic  temperatures  because  the 
recorded  information  is  actually  stored  as 
population  modulation  of  hyperfine  levels 
in  the  ground  state  (6).  Because  the  grating 
is  a  direct  result  of  resonant  excitation  and 
requires  no  secondary  process,  its  formation 
is  instantaneous  and  requires  modest  laser 
power;  the  data  recording  is  therefore  in¬ 
herently  fast  and  efficient,  giving  rise  to 
high  data  transfer  rates,  fast  access  times, 
and  short  latency.  These  features  are  ideal 
for  applications  such  as  a  dynamic  memory. 

Data  stored  as  population  gratings  are 
retrieved  by  illuminating  the  sample  with 
another  laser  pulse  known  as  the  read  puke 
(Fig.  2),  which  reconstructs  the  object  field 
as  in  conventional  holography.  However, 
unlike  in  conventional  holography,  the  re¬ 
constructed  field  appears  as  an  emitted  laser 
puke  and  is  temporally  distinct  from  the 
read  pulse,  just  as  in  case  of  recording.  For  a 
given  temporal  separation  t  between  the 
write  and  data  pulses,  the  reconstructed 
pulse  will  occur  a  time  t  after  the  arrival  of 
the  read  puke  (Fig.  1)  and  has  spatial  fea¬ 
tures  identical  to  those  of  the  stored  data 
pulse.  This  delayed  radiation  is  often  re¬ 
ferred  to  as  photon  echoes.  Mathematically, 
the  Fourier  transforms  of  the  echo  and  three 
input  puke  electrical  fields  £  are  related  in 
the  following  manner 

£e(v)  *  £^(v)£d(v)£r(v)  (l) 

where  the  subscripts  e,  W,  D,  and  R  denote 
echo,  write,  data,  and  read  pukes,  respec¬ 


tively,  and  v  denotes  the  frequency.  Be¬ 
cause  the  echo  is  a  replica  of  the  data  puke 
both  spatially  and  temporally,  the  peak 
readout  rate  (the  number  of  bits  per  echo)  is 
identical  to  the  peak  recording  rate  (the 
number  of  bits  per  data  pulse)  in  a  time- 
domain  digital  memory. 

The  material  we  chose  to  demonstrate  this 
recording  scheme  was  crystalline  Eu3+: 
Y2Si05  (0.1  atomic  %,  7  mm  thick).  We  used 
its  %-5Do  transition  at  site  1  ( — -579.88  nm), 
which  has  an  inhomogeneously  broadened 
linewidth  of  ~4  GHz  and  a  dephasing  time  of 
—2  ms  (9).  The  product  of  the  two  parameters 
gives  the  theoretical  upper  bound  for  the  stor¬ 
age  capacity  per  spatial  location  as  a  result  of 
the  presence  of  the  wavelength  dimension. 
This  material  has  two  optically  equivalent 
sites  for  the  transition;  the  other  site  (site  2) 
has  a  comparable  linewidth  and  dephasing 
time  (9),  making  the  material  a  good  candi¬ 
date  for  high-capacity  applications.  During 
the  experiments,  the  crystal  was  kept  at  a 
temperature  of  about  4  K  in  a  flowing  He- 
vapor  cryostat. 

Digital  data  pages,  each  of  which  con¬ 
tained  3360  bits  of  data,  were  generated  by 
illuminating  a  spatial  light  modulator 
(SLM)  (10)  with  laser  pukes  as  short  as 


11.2  |xs  (II),  giving  rise  to  a  maximum  peak 
(burst)  recording  rate  of *300  megabits  per 
second.  Each  page  was  stored  in  a  distinct, 
narrow  ( — 500  kHz)  frequency  channel  (12) 
within  the  inhomogeneously  broadened  ab¬ 
sorption  line  by  wavelength-division  multi¬ 
plexing  (WDM).  The  separation  between 
two  adjacent  channels  was  chosen  to  be ' 
'■'■*800  kHz.  There  are  distinct  advantages  to 
WDM  (13):  (i)  Recording  and  retrieval  can 
be  accomplished  without  the  use  of  any 
moving  parts  for  laser  beam  steering,  so  that 
it  is  inherently  reliable  and  often  provides 
short  access  times  (12);  (ii)  there  is  no 
significant  effect  on  diffraction  efficiency; 
and  (iii)  each  data  page  can  be  recorded  in 
an  identical  manner,  eliminating  the  need 
for  scheduled  or  incremental  recording,  as  is 
used  in  angular-division  multiplexing  (14). 

The  retrieved  data  were  captured  by  an 
intensified  charge-coupled  device  (CCD) 
camera  and  subsequently  digitized  by  a  frame 
grabber  (15).  The  SLM  and  CCD  operated  at 
the  video  rate  (30  frames  per  second,  which  is 
the  maximum  frame  rate  allowed  for  the  de¬ 
vices),  and  the  resulting  digitized  data  were 
transferred  at  the  same  speed  to  a  computer’s 
dynamic  memory.  At  the  present  time,  the 
SLM  limits  the  sustainable  recording  speed  to 
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Fig.  2.  Schematic  of  the 
mechanism  for  the  formation 
of  the  population  grating  and 
echoes  under  the  laser  exci¬ 
tations.  Symbols  Ig)  and  le> 
denote  the  ground  and  excit¬ 
ed  states  of  an  electronic 
transition,  respectively,  and 
the  vertical  labeled  f  repre¬ 
sents  the  frequency. 
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-lie  kilobits  per  second,  substantially  lower 
than  the  burst  rate  of  300  megabits  per  second 
permitted  bv  the  storage  material;  however, 
this  limit  is  not  fundamental.  Hiizh-speed 


Fig.  3.  (A)  Input  data  image  after  transmission 
through  the  Eu3*:Y2Si05  sample.  (B)  Recon¬ 
structed  image  after  storage  in  Eu3~:Y2Si05  for 
200  pis.  (C)  Image  reconstructed  under  the  same 
conditions  as  in  (B)  but  with  weak  overall  intensity. 
(D)  Image  reconstructed  5  min  after  the  comple¬ 
tion  of  recording.  In  (A)  through  (C),  the  length  of 
the  data  puise  was  1 1 .6  fis.  and  that  in  (D)  was  50 
\ls.  These  times  correspond  to  I/O  rates  of  300 
and  67  megabits  oer  second,  respectively. 
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SLMs  with  a  full  frame  rare  of  3  kH:  are 
commercially  available;  ones  with  higher 
speeds  are  being  considered  for  development. 
Semiconductor  laser  arrays  may  be  used  in  the 
future  to  compose  binary  digital  pages  without 
the  need  of  SLMs.  Such  devices  will  have  a 
single-channel  modulation  rate  sufficiently 
high  (in  excess  of  1  GH:)  to  match  the  burst 
rate  sec  by  the  page  exposure  time. 

In  digital  memories,  one  important  param¬ 
eter  determining  the  memory’s  performance  is 
the  bit  error  rate  (BHR),  which  measures  how 
often  an  error  bit  occurs  on  average.  This 
parameter  is  closely  related  to  other  perfor¬ 
mance  parameters  of  the  system,  such  as  stor¬ 
age  capacity,  data  transfer  rates,  and  data  ac¬ 
cess  times.  In  general,  for  a  given  system, 
attempting  to  reduce  the  BER  leads  to  a  de¬ 
crease  of  storage  capacity  or  data  transfer  rate, 
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Fig.  4.  (A)  Histogram  showing  the  intensity  distri¬ 
butions  of  0’s  and  l ’s  that  were  retrieved  from  the 
time-domain  memory.  (B)  Normalized  histogram 
obtained  after  applying  the  local  threshold  tech¬ 
nique.  In  obtaining  these  results,  we  used  a  peak 
recording  rate  of  300  megabits  per  second.  The 
100  pages  were  wavelength -division  multiplexed 
at  a  single  storage  cell  measuring  about  1 .0  mm 
by  1 .0  mm  by  7.0  mm  at  a  rate  of  30  frames  per 
second  with  a  total  recording  time  of  -3.3  s.  All  of 
the  pages  used  here  were  obtained  from  the  same 
experiment.  For  convenience,  each  page  used  for 
BER  analysis  was  retrieved  after  being  transferred 
to  the  memory  for  200  ^s.  The  use  of  short  stor¬ 
age  time  prevented  the  undesired  page-to-page 
fluctuation  caused  by  the  slow  frequency  drift  of 
our  laser.  An  example  of  long-term  storage  is 
shown  in  Rg.  3D.  The  echo  (or  “diffraction")  effi¬ 
ciency  for  the  short-term  storage  was  measured 
to  be  -5  x  10-3. 


or  both,  because  a  low  BER  requires  high 
signal-to-noise  ratios,  which  implies  the  use  of 
large  bit  dimensions  in  the  case  of  parallel 
recording.  For  a  given  frame  size,  it  translates 
into  a  decrease  of  the  data  transfer  rate. 

For  an  optical  storage  system  operating 
at  high  I/O  speeds  with  low  laser  power 
illumination,  the  retrieved  data  are  often  at 
a  signal  level  that  is  limited  by  shot  noise, 
which  is  a  result  of  the  quantum  nature  of 
photons.  For  page-formatted  recording,  the 
process  is  further  complicated  by  the  non- 
uniform  illumination  of  an  SLM.  Because  of 
these  undesired  conditions  associated  with 
a  low  photon  budget,  such  a  svstem  often 
exhibits  large  page-ro-page  fluctuations  and 
intrapage  variations. 

Figure  3A  shows  a  sample  input  page 
used  in  data  recording.  Spatial  variation  of 
the  intensity  is  apparent.  This  variation 
reflects  the  spatial  beam  mode  of  the  ad¬ 
dressing  laser,  contributing  partially  to  the 
spatial  signal  variations  within  a  page  dis¬ 
cussed  earlier.  The  ocher  contributions  in¬ 
clude  shot  noise. 

The  page-to-page  fluctuation  is  often 
more-  noticeable  in  the  retrieved  data  be¬ 
cause  it  is  magnified  by  less-than-ideal  re¬ 
cording  and  readout  conditions.  Figure  3,  B 
and  C,  show  rwo  retrieved  pages  that  rep¬ 
resent  the  brightest  and  dimmest  among  a 
total  of  100  pages  recorded  in  an  experi¬ 
ment.  They  are  single  readout  events  with 
no  signal  averaging.  The  intensity  differ¬ 
ence  between  the  two  pages  is  >300%, 
which  is  mainly  caused  by  the  frequency 
instability  of  the  laser  system  (12). 

Examination  of  the  histogram  of  all  100 
pages  (totaling  336  kilobits  of  data)  reveals 
that  the  intensity  distributions  of  the  l’s  (or 
ON’s)  and  0’s  (or  OFF’s)  partially  overlap 
(Fig.  4A).  If  a  single  threshold  were  used  to 
determine  all  of  the  data  bits,  it  would  lead 
to  a  BER  as  large  as  10“ \  which  is  not 
suitable  for  any  storage  application. 

The  broad  distribution  curves  in  Fig.  4 A 
are  mainly  a  reflection  of  the  interpage  inten¬ 
sity  fluctuations  and  intrapage  variations.  The 
broadening  also  includes  cross  talk  between 
neighboring  l’s  and  0’s.  If  these  intensity  vari¬ 
ations  and  cross  talk  were  removed,  the  dis¬ 
tributions  for  both  l’s  and  0’s  would  follow 
closely  scaled  Poisson  distributions,  which  has 
been  confirmed  experimentally  for  known  im¬ 
ages.  Under  this  condition,  a  single  threshold 
can  be  obtained  for  the  whole  page  with 
knowledge  of  the  average  of  the  intensities  of 
l’s  and  0’s  only.  This  finding  allows  us  to 
develop  the  following  scheme  to  rapidly  and 
accurately  read  page-formatted  binary  data 
that  have  large  intrapage  intensity  variations 
and  large  interpage  fluctuations. 

We  use  recently  determined  data  bits  to 
obtain  a  local  threshold  for  the  next  bit 
under  consideration.  The  local  region  is 
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chosen  in  such  a  way  that  it  occupies  an 
area  that  is  small  enough  so  that  the  illu¬ 
mination  is  approximately  uniform  yet  large 
enough  to  have  a  sufficient  number  of  bits 
for  an  adequate  average.  Under  such  a  con¬ 
dition,  a  threshold  can  be  derived  from  the 
scaled  Poisson  distributions  to  accurately 
determine  the  state  of  the  bit,  provided  that 
the  intensity  within  a  page  varies  smoothly. 
We  adopted  the  following  formula  for  the 
threshold 

Jthr  “  f 0  P  \  fo 

+  4Ui -PnM  Mfo  +  Pv^)]  (2) 

where  I0  and  fj  are  the  average  intensities  of 
O’s  and  Ts,  q  is  a  weighting  factor  (chosen 
empirically  to  be  l/3),  and  p  is  a  measure  of 
at  how  many  standard  deviations  above  I0 
and  below  we  can  place  the  threshold. 
Here,  we  use  p  =  11.  The  optimum  num¬ 
bers  used  in  averaging  were  found  to  be  4 
for  O’s  and  8  for  l’s. 

Initial  values  of  l0  and  are  needed  to 
jump  start  the  analysis  of  each  page.  This  is 
accomplished  by  permanently  setting  to  0 
four  known  bits  at  the  comers  of  a  4  X  4 
region  in  the  center  of  the  SLM,  where  the 
signal- to-noise  rate  is  largest.  We  average 
the  4  bits  to  estimate  J0.  A  crude  estimate  for 
1 1  can  be  obtained  by  averaging  the  intensi¬ 
ties  of  the  remaining  12  unknown  bits.  The 
threshold  for  the  block  is  then  calculated 
from  Eq.  2  by  assuming  the  block  intensity 
Iblock  =  »o  +  fj)/2»  which  implies  an  equal 
number  of  l’s  and  0's  in  the  block.  This 
approach  yields  a  BER  of  10“6  for  the  12-bit 
central  region,  or  equivalently,  a  10“ 3  prob¬ 
ability  of  having  an  error  bit  per  page. 

We  then  continue  this  analysis  process  by 
spiraling  outward  from  the  center  of  a  page 
toward  the  edges,  with  the  intensity  decreas¬ 
ing  gradually  as  the  analysis  proceeds.  The 
neighboring  l’s  and  O’s  are  stored  in  two 
separate  register  arrays  and  updated  in  a  first- 
in-first-out  manner  to  obtain  their  running 
averages.  Thus,  each  data  bit  is  determined  in 
real  time  with  the  use  of  Eq.  2.  In  comparison 
to  existing  schemes  for  reading  page-format¬ 
ted  digital  data  (2,  16),  this  automatic  thresh¬ 
old  calibration  technique  uses  no  coding  bits 
and  requires  no  mathematical  iteration.  These 
advantages  lead  to  a  significant  increase  in 
the  readout  rate,  which  is  ideal  for  high-speed 
applications. 

Using  this  approach,  we  repeated  the 
above  experiment  and  obtained  a  histogram 
for  the  same  100  pages  used  earlier  (Fig.  4B). 
This  histogram  was  normalized  to  the  respec¬ 
tive  threshold  of  each  data  bit  and  shows  two 
clearly  separate  distribution  curves.  The  sep¬ 
aration  is  —40%  of  the  full  width  of  the 
distribution  for  O’s.  This  large  separation  sug¬ 
gests  an  effective  increase  of  the  system’s  dy¬ 
namic  range  by  -“40%.  The  raw  BER  is  esti¬ 


mated  to  be  —10“'.  We  further  verified  this 
number  by  repeating  the  experiment  10  times: 
The  two  distribution  curves  in  the  resulting 
histogram  remained  separate. 

According  to  our  BER  measurement,  our 
current  storage  system  is  capable  of  recording 
close  to  10  megabits  of  digital  data  without 
generating  an  error  bit.  To  demonstrate  it,  we 
stored  a  high-resolution  black-and-white  pho¬ 
tograph  of  Albert  Einstein  (Fig.  5A).  The 
digitized  photograph  contains  106  pages,  to¬ 
taling  356  kilobits,  of  digital  data.  Figure  5B 
shows  the  reconstructed  photograph  after  be¬ 
ing  copied  to  and  read  back  from  the  memory. 
The  enure  photograph  was  stored  at  a  single 
spatial  location  in  Eu3TY:Si05  using  WDM 
with  a  peak  recording  speed  of  300  megabits 
per  second,  and  no  error  correcting  code  of 
any  kind  was  used.  The  total  material  band¬ 
width  used  was  only  —85  MHz. 

Because  Eu3TY2Si05  has  two  equiva¬ 
lent  optical  sites  with  a  total  of  —8  GHz 
inhomogeneous  bandwidth,  85  MHz  repre¬ 
sents  only  —1%  of  the  available  storage 
capacity  at  this  location.  In  principle,  100 
such  photographs  can  be  wavelength  mul¬ 
tiplexed  at  this  spatial  location,  which  im¬ 
plies  a  total  raw  capacity  of  —4.5  mega¬ 
bytes.  Taking  into  account  the  active  stor¬ 
age  volume  used,  the  projected  density  is 
—6  X  108  bits  per  cubic  centimeter. 

The  above  estimates  were  based  on 
simple  extrapolation  of  our  current  exper¬ 
imental  results,  which  were  obtained  un¬ 
der  conditions  that  were  far  from  optimum 
(15);  therefore,  these  numbers  do  not  rep¬ 
resent  the  fundamental  limits  of  the  stor- 
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age  system.  The  projected  capacity  can  be 
increased  substantially  by  using,  for  exam¬ 
ple,  a  matched  SLM,  CCD,  and  frame 
grabber  (3).  Furthermore,  the  use  ot  SLMs 
with  lower  insertion  loss  would  also  in¬ 
crease  the  storage  density  because  higher 
spatial  resolution  can  be  achieved  with  a 
high-power  addressing  beam.  Currently, 
SLMs  with  an  insertion  loss  of  <20%  are 
commercially  available.  With  the  rapid 
advance  of  peripheral  devices,  such  as 
SLMs  and  CCD  cameras,  the  achievable 
storage  capacity  is  likely  to  be  several 
orders  of  magnitude  greater  than  what  is 
projected  here. 

The  recording  and  readout  speeds  would 
also  be  increased  by  the  use  of  a  SLM  with 
lower  insertion  loss.  We  estimate  that  the 
length  of  the  data  pulse  can  be  reduced  by  a 
factor  of  5  without  loss  of  signal  strength  by 
using  a  SLM  with  20%  insertion  loss.  This 
adjustment  would  allow  an  I/O  rate  of  —  1 .5 
gigabits  per  second.  Furthermore,  by  using 
matched  SLMs  and  CCDs,  an  I/O  rate  in 
excess  of  10  gigabits  per  second  should  be 
achievable  without  any  penalty  in  raw  BER. 

Parallel  digital  data  storage  making  use  of 
the  time-domain  approach  could  be  devel¬ 
oped  into  a  high-performance  page-oriented 
dynamic  optical  memory.  As  the  technology 
for  parallel  data  access  becomes  more  mature, 
this  time-domain  holographic  approach  along 
with  other  optical  holographic  storage  tech¬ 
niques  (such  as  photorefiractive  memory)  may 
fulfill  their  long-awaited  promise  for  high¬ 
speed,  high-density  data  storage.  A  number  of 
issues  still  need  to  be  addressed  before  this 


Fig.  5.  (A)  Digitized  photograph  displayed  using  its  original  binary  data.  (B)  Same  photograph  obtained 
after  the  binary  data  was  transferred  to  and  then  recalled  back  from  the  time-domain  memory.  [Reprint¬ 
ed  with  permission  from  UPl/Corbis-Bettman] 
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